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The 


Preface 


purpose  of  this  thesis  effort  was  to  find  the 
resonant  frequencies  of  the  Compressor  Research  Facility 
located  at  Wright-Patterson  AFB,  Ohio.  The  general  approach 
to  the  study  was  to  treat  the  facility  as  a  fluid  transmission 
line  and  solve  for  the  transfer  function  of  the  facility.  The 
solution  is  analogous  to  electrical  theory  where  the  gain  and 
phase  shift  of  an  alternating  input  signal  is  determined  at 
the  output  end  of  the  line. 

I  had  the  good  fortune  to  be  advised  by  Dr.  Milton 
Franke,  who  has  established  an  outstanding  reputation  for  his 
work  in  fluid  transmission  lines,  as  well  as  in  other  areas  of 
mechanical  and  aeronautical  engineering.  Dr.  Franke 's 
guidance  through  the  project  was  the  primary  reason  this 
thesis  effort  was  successful. 

This  project  involved  both  theoretical  and  experimental 
efforts  and  the  latter  could  not  have  been  accomplished 
without  the  help  of  Andy  Pitts  of  the  AFIT  Aeronautical 
Engineering  Laboratories.  Andy's  help  with  setting  up  the 
transducers,  plumbing,  data  acquisition  system,  and  pneumatic 
valve  was  essential  to  the  success  of  the  experiment.  Also, 
Mr.  Jan  LeValley  and  Mr.  David  Driscoll  at  the  AFIT  model  shop 
constructed  the  CRF  model  with  superb  accuracy,  resulting  in 
excellent  experimental  data. 
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Research  Facility,  was  greatly  appreciated.  Although  often 
busy  with  the  daily  operation  of  the  CRF,  Doug  always  took  the 
time  to  answer  any  questions  I  had  concerning  the  facility. 

Also  greatly  appreciated  is  the  support  and  team-oriented 
efforts  of  Capt  Brian  Hull.  Though  busy  with  his  own  thesis, 
Brian  worked  with  me  to  connect  his  control  system  design  to 
my  theoretical  output.  His  faith  in  nty  work  was  crucial  to  my 
self-esteem  as  well  as  to  the  usefulness  of  the  results  in 
this  thesis. 

I  would  like  to  thank  Maj  Tom  Buter  and  Capt  John  Doty, 
readers  of  this  thesis,  for  their  support  of  the  project  and 
helpful  suggestions  along  the  way.  Their  effort  to  understand 
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to  other  engineers . 
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Abstract 


The  response  of  the  Wright  Laboratory  Compressor  Research 
Facility  to  a  small  amplitude,  acoustic,  sinusoidal 
disturbance  was  investigated.  The  fluid  transmission  line 
equations  for  laminar,  one-dimensional  flow  in  a  circular 
duct  were  s<^Tved  and  verified  through  a  laboratory 
experiment  using  a  scale  model  of  the  facility.  The 
resonant  frequencies  of  the  facility  were  determined  for  a 
variety  of  flow  conditions.  Techniques  for  analytically 
modeling  end  impedance  and  flow  straightening  elements  were 
also  studied.  The  fundamental  frequency  of  the  facility  was 
determined  to  be  between  5.5  and  6.5  Hz  depending  on  the 
flow  conditions  and  geometry  configurations  specified. 
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DYNAMIC  RESPONSE 


OF 

A  COMPRESSOR  RESEARCH  FACILITY 

I.  Introduction 

The  Study  of  the  dynamic  response  of  the  Air  Force  Wright 
Laboratory  Compressor  Research  Facility  (CRF)  was  proposed  for 
the  purpose  of  determining  the  resonant  frequencies  of  the 
facility.  Currently,  a  new  inlet  valve  control  system  design 
is  under  development  in  a  parallel  effort  (Hull: 1993).  This 
control  system  design  requires  accurate  resonant  frequency 
input  in  order  to  perform  properly.  This  study  determines  the 
resonant  frequencies  through  a  modeling  of  the  CRF  as  a  fluid 
transmission  line. 

Fluid  transmission  lines  have  been  studied  in  depth  for 
over  40  years.  In  the  early  years,  these  studies  were 
conducted  to  yield  baseline  equations  and  data  for  design  of 
fluidic  components  to  be  used  as  amplifiers,  diodes,  switches, 
and  other  circuit  elements.  The  early  studies  resulted  in  a 
sound  knowledge  base  of  fluid  transmission  line  theory 
( Iberall : 1950 ;  Nichols : 1962 ;  Brown:1962;  Kirshner  and 

Katz: 1975).  In  addition,  efforts  to  gather  baseline  data 
resulted  in  theory  verification  for  many  types  of  fluid 
transmission  lines  (Krishnayer  and  Lechner :1967 ;  Bergh  and 
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Tijdeman : 1965 ;  Franke:1969,  1979,  1972) .  The  theory  and  some 
experimental  data  from  the  early  work  apply  to  the  solution  of 
the  CRF  dynamic  response. 

1.1.1  Fluid  Transmission  Lines.  The  theory  centers  on 
the  idea  that,  for  a  small  amplitude  input  signal,  fluid 
transmission  lines  are  fundamentally  similar  to  electrical 
transmission  lines.  To  begin  with,  the  fluid  line  is 
comprised  of  the  same  passive  elements  in  an  electrical  line; 
namely,  resistance,  conductance,  inductance,  and  capacitance. 
Viscous  dissipation  causes  losses  in  the  fluid  line  in  the 
same  manner  as  resistance  and  conductance  cause  losses  in 
potential  in  an  electrical  line.  Conpressed  fluid  stores 
energy  in  the  same  manner  as  a  capacitor  in  an  electrical 
line.  Boundary  layer  interactions  cause  a  phase  shift  in  the 
signal  in  the  same  manner  as  an  inductor  in  an  electrical 
line.  Further,  the  pressure  difference  between  the  exit  and 
inlet  of  a  fluid  line  is  analogous  to  voltage  while  the 
volumetric  flow  is  considered  the  fluid  equivalent  of  current. 
A  circuit  made  up  of  passive  elements  has  resonant  frequencies 
and  at  these  frequencies  the  output  can  be  much  larger  than 
the  input,  resulting  in  a  large  gain.  Given  a  sinusoidal 
input  signal,  the  gain  and  phase  shift  of  the  signal  at  the 
output  of  the  circuit  can  be  calculated  if  the  resistance, 
conductance,  capacitance,  and  inductance  are  known.  For  a 
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sinple  electrical  circuit  design  in  which  the  passive  elements 
are  specified,  finding  the  resonant  frequencies,  the 
associated  gain,  and  corresponding  phase  shift  is  quite 
tractable.  For  transmission  lines;  however,  considerable 
effort  is  required  as  values  for  resistance,  conductance, 
inductance,  and  capacitance  must  be  theoretically  determined. 
Fortunately,  these  parameters  can  be  found  for  the  fluid  line 
by  solving  equations  developed  by  Iberall  (1950:85-108)  and 
Nichols  (1962:5-14). 

In  this  study,  Nichols'  work  is  used  to  determine  the 
transmission  line  characteristics  of  the  CRF.  In  particular, 
the  study  focuses  on  finding  the  resonant  frequencies  of  the 
facility. 

1.1.2  The  Coi  jressor  Research  Facility.  The  CRF  is  a 
unique  facility  used  to  test  experimental  jet  engine 
compressors.  Flow  is  introduced  through  five  large  butterfly 
valves  at  the  inlet,  drawn  through  several  flow  conditioning 
elements,  and  exhausted  through  the  test  article  (compressor)  , 
The  flow  conditions  are  varied  according  to  a  test  plan 
developed  for  each  compressor  that  typically  calls  for  several 
runs  at  different  flow  conditions.  Often,  specified  flow 
conditions  are  met  by  setting  a  desired  static  pressure  at  the 
face  of  the  compressor;  thereby  simulating  a  predetermined 
altitude  for  testing. 
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Operators  of  the  CRF  primarily  control  the  pressure  at 
the  compressor  inlet  through  the  manipulation  of  the  five 
butterfly  valves.  However,  when  the  speed  of  the  test  article 
is  changed,  the  static  pressure  at  the  conpressor  face  also 
changes.  Due  to  this  effect,  it  is  not  a  simple  process  to 
evaluate  and  test  a  compressor  at  different  speeds  and 
constant  pressure.  Presently,  a  feedback  control  system  is 
being  designed  to  read  the  static  pressure  change  caused  by  a 
change  in  the  compressor  speed  and  accordingly  adjust  one  or 
more  of  the  five  inlet  valves  to  hold  a  constant  pressure  at 
the  inlet  of  the  test  article. 

However,  when  the  control  system  adjusts  the  valves,  the 
steady  inlet  flow  will  be  disturbed.  As  a  result,  unsteady 
pressure  waves,  at  a  frequency  corresponding  the  valve 
movement  rate,  will  be  produced  and  transmitted  down  the 
facility.  These  waves,  should  they  be  generated  at  one  of  the 
resonant  frequencies  of  the  facility,  could  become  quite  large 
and  ultimately  destroy  the  test  article. 

The  control  system,  then,  must  be  designed  to  avoid 
operation  at  the  fundamental  frequencies.  The  purpose  of  this 
study  is  to  prc/j.de  that  frequency  information  to  the  controls 
designer  so  the  control  system  can  be  developed  with 
confidence  that  the  resonant  frequencies  will  be  avoided. 
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The  primary  objective  of  this  study  is  to  determine  the 


resonant  frequencies  of  the  CRF.  The  frequencies  are 
determined  from  the  fluid  transmission  line  theory  presented 
by  Nichols  (1962).  The  Nichols  equations  are  modified  to 
account  for  the  specific  geometry  of  the  facility  and  then 
solved  for  various  mass  flows  and  compressor  radii.  The 
solution  process  involves  meeting  the  following  specific 
objectives : 

1.  Write  and  verify  a  computer  program  that  determines 
the  frequency  response  of  the  CRF; 

2.  Design  and  test  a  scale  model  of  the  CRF  and  verify 
the  analytical  results; 

3.  Determine  the  correct  analytical  model  of  the  screens 
and  honeycomb  elements  in  the  CRF;  and 

4.  Obtain  the  transient  response  of  the  laboratory 
model . 

1.3  Thesis  Outline 

The  entire  discussion  is  comprised  of  eight  chapters. 
First,  the  governing  equations  for  fluid  transmission  lines 
are  detailed  in  Chapter  II.  The  physical  systems  are  outlined 
in  Chapter  III,  including  the  various  configurations  of  the 
CRF  and  the  laboratory  model .  The  experimental  apparatus  and 
data  acquisition  equipment  are  described  in  Chapter  IV. 


Chapter  V  outlines  the  computer  program  written  to  solve  for 
the  frequency  response  of  the  CRF  and  laboratory  model .  The 
results  of  the  program,  as  applied  to  the  CRF  are  discussed 
and  presented  in  Chapter  VI.  The  experimental  results  from 
the  scale  model  are  presented  in  Chapter  VII. 

Conclusions  on  these  results  are  discussed  in  Chapter 
VIII  along  with  recommendations  for  future  studies  of  this 
type.  Finally,  the  results  of  several  intermediate  studies 
are  presented  in  Appendices . 
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In  this  chapter,  the  basics  of  pipe  resonance  are 
covered,  followed  by  the  circuit  theory  for  fluid  transmission 
lines,  and  a  summary  of  the  governing  equations.  The 
transmission  line  equations  for  a  constant  diameter  circular 
line  with  a  small,  acoustic  input  signal  are  presented  and 
modifications  to  include  velocity  effects  discussed.  The 
chapter  is  concluded  with  equations  for  more  complex 
geometries,  including  cascaded  and  branched  lines. 

2.2  Basic  Pipe  Resonance  Equations 

The  fluid  transmission  line  theory  is  really  an  extension 
of  the  basic  equations  for  a  resonating  pipe,  such  as  an  organ 
pipe.  For  a  pipe  of  length,  1,  and  air  with  speed  of  sound, 
a,  the  resonant  frequencies  are  given  by  (Haliday  and  Resnick, 
1981:326)  : 

f  =  .££ 

r«aa,b  (2-1) 

n  —  1,3,5,... 

for  a  blocked  pipe.  For  an  open  pipe,  the  resonant 
frequencies  are: 
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(2-2) 


f  _  na 
-  Yl 

n  =  1,2,3,... 

From  Equations  (2-1)  and  (2-2),  the  following 
observations  are  made.  First,  the  resonant  frequencies  are 
inversely  proportional  to  the  length  of  the  pipe.  As  such, 
long  pipes  will  have  lower  resonant  frequencies  than  short 
pipes.  Second,  for  two  pipes  of  the  same  length,  one  open  and 
one  blocked,  the  blocked  pipe  resonant  frequencies  will  be  one 
half  of  the  open  pipe  resonant  frequencies.  Therefore,  for 
the  same  length,  the  open  pipe  will  have  higher  resonant 
frequencies  than  the  closed  pipe.  Third,  depending  on  the 
frequency  range,  length  of  pipe,  and  speed  of  sound,  the 
blocked  pipe  can  have  more  resonant  frequencies  in  the  range 
than  the  open  pipe.  The  latter  notion  is  more  fully  developed 
in  Chapter  III. 

2.3  Circuit  Theory  for  Fluid  Transmission  Lines 

The  circuit  representation  for  a  typical  fluid 
transmission  line  with  passive  elements  and  a  small,  acoustic 
input  signal  is  illustrated  in  Figure  2.1.  Pj  is  the  input 
pressure  signal  and  Q  is  the  flow  rate.  At  some  point  along 
the  line,  the  pressure  signal  is  detected  as  Pj  and  its  value 
depends  on  the  load  impedance  Z^,  as  well  as  on  the  line 
impedance,  Z,  and  admittance,  Y. 
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For  any  given  circular  fluid  line,  knowledge  of  Z,  Y,  and 
Zl  allows  solution  of  the  system  gain,  Pj/Pj,  and  phase  shift, 
0,  at  any  frequency.  The  result  fully  defines  the  dynamic 
response  of  the  transmission  line  for  a  small  sinusoidal  input 
signal . 

2.4  Governing  Equations 

Analytical  methods  for  determining  Z,  Y,  and  Pi/Pi  for  the 
fluid  transmission  line  have  been  developed  from  the  Navier- 
Stokes  equations  (Nichols,  1962:5-14).  The  solution  of  the 
Navier-Stokes  equations  is  reduced  in  conplexity  according  to 
the  following  simplifying  assunptions : 

1.  one-dimencional  plane  wave; 

2.  axial  symmetry; 

3.  ideal  gas; 

4.  small  acoustic  signals; 

5.  laminar  flow; 


Figure  2.1.  Fluid  Transmission  Line  Circuit  Diagram 
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6.  first  order  approximations; 

7.  negligible  axial  velocity  and  tenperature  gradients; 

8.  uniform  pressure  at  each  cross  section;  and 

9.  no  throughflow. 

The  simplified  momentum,  energy,  and  continuity  equations  are 
(Kirshner  and  Katz,  1975:76): 

^  {momentum)  <2-3) 

^  p  ^  r  or\ 


PC.|I 


(energy)  (2-4) 


+  o  = 


(continuity)  (2-5) 


Nichols  (1962:9)  has  shown  these  equations  reduce  to  the  wave 
equation  if,  at  the  inlet  (x  =  0),  the  alternating  pressure 
and  alternating  flow  rate  are  specified: 


=  ZYP 


(2-6) 


In  the  wave  equation,  the  impedance,  Z,  and  the  admittance,  Y, 
are  given  by 


2  Jj  (8(0/0^,) 

j3/2  (8C0/(a,,)  1/2] 


(2-7) 


(2-8) 


1  +  2(Y  -  1) 

Where  the  characteristic  frequency  for  a  circular  line  is 
given  by: 


(2-9) 


The  three  parameters  that  completely  characterize  the 
transmission  line  are  the  end  iitpedance,  Z^,,  the 
characteristic  impedance,  Zc,  and  the  propagation  operator,  F. 
The  end  impedance  is  conprised  of  a  resistance  and  inertance 
and  is  discussed  in  detail  in  the  next  section.  Z^  and  F  are 
dependent  on  Z  and  Y  as 

»  iZ/Y)^'^  (2-10) 

F=(ZY)i'2  .  (2-11) 

Finally,  Kirshner  &  Katz  (1975:68)  give  the  gain  of  the  line 
as : 


Pi  ^  {Z^/ZJ 

Pj  sinh(FJ)  +(Zj^/Z^)cosh(rj) 


(2-12) 


and  the  phase  angle  as: 
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0  =  arctan 


Iin(Pj/P^) 

Re(Pj/Pi) 


(2-13) 


Given  the  flow  conditions  and  geometry  of  the  fluid  line,  the 
only  remaining  variable  to  be  specified  is  the  end  impedance, 
Zl- 


2.4.1  End  Impedance .  The  end  impedance  is  comprised  of 
a  resistance,  R,  and  an  inertance,  L: 

=  R  +  j(OL  .  (2-14) 


In  the  case  of  a  blocked  line,  is  infinitely  large.  On  the 
other  hand,  for  an  open  line,  is  zero.  For  an  orifice- 
terminated  line,  Zl  may  be  represented  as  (Kirshner  and  Katz, 
1975:16)  : 


(2-15) 


where  A  is  the  area  of  the  smallest  cross  section  and  Kl  is 
the  loss  coefficient  of  the  orifice. 

According  to  Kirshner  &  Katz  (1975:56),  inertance  is 
unimportant  unless  the  orifice  length  is  short  in  comparison 
to  the  wavelength  of  the  input  signal .  When  required,  the 
inertance  of  a  line  is  expressed  as: 


(2-16) 
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where  1  and  A  are  the  length  and  area  of  the  orifice, 
respectively . 

For  a  single,  constant  diameter  line  with  no  throughflow, 
the  above  equations  satisfactorily  describe  the  dynamic 
response  of  a  system.  However,  the  end  impedance  for  an 
orifice-terminated  line  is  strongly  dependent  on  mass  flow 
rate.  As  such,  it  is  natural  to  extend  the  expressions  for 
the  propagation  operator  and  characteristic  impedance  to 
include  the  effects  of  throughflow  as  well. 

2.4.2  Throughflow  Effects.  If  throughflow  is  present,  an 
adjustment  can  be  made  to  the  characteristic  impedance  and  the 
propagation  operator.  Katz,  Housner,  and  Eisenberg 
(1974:278-281)  derived  approximations  to  F  and  Z^.  that 
incorporate  velocity  effects. 

With  these  approximations,  the  propagation  operator  is 
adjusted  for  throughflow  as  follows: 


r. 


r  (1  + 

1  - 


(2-17) 


where 


m 


M 

2  Z, 


z 


ca 
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The  gain  is  then  calculated  as: 


Pi  ,  e'-Ml  + 

Pi  “  (1  +  jn2)  1/2 cosh (Fyl)  +  [(Z^/ZJ  -  jn]  sinh(r^I) 


(2-18) 


with 


b  =  "Xg , 

2(1-  M') 


The  characteristic  impedance  then  becomes: 


=  Z^[(l  +  in2)i/2  -  n,] 


(2-19) 


2.4.3  Cascaded  Lines.  The  theory  for  a  single,  constant 
diameter  line  is  easily  extended  to  a  system  of  cascaded  lines 
with  different  diameters.  In  a  system  of  cascaded  lines,  the 
overall  gain  is  found  by  multiplying  the  gains  of  the 
individual  lines.  Further,  a  new  variable  is  introduced,  Zj, 
as  the  input  impedance  of  a  line.  In  general,  Zi  is  equal  to 
the  end  impedance  of  the  upstream  line  and  is  computed  through 
(Kirsliner  and  Katz,  1975:67)  : 


= 


Z^cosh(ri)  +Z^sinh(ri) 
Z^sinh(rj)  +Z^cosh(ri) 


(2-20) 
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Often,  cascaded  line  behavior  is  significantly  different 
from  the  behavior  of  a  single  diameter  line  of  the  same 
length.  Depending  on  the  length  of  the  lines  in  the  cascade, 
waves  reflecting  at  the  boundaries  can  change  the  resonant 
frequencies  from  those  of  a  constant  diameter  line.  Thus,  the 
overall  transfer  function  can  change  when  the  cascade  is 
introduced. 


2.4.4  Branched  lines.  A  branched  line  allows  for  a 
second  flow  path  and  can  also  significantly  alter  the  transfer 
function  of  a  system.  Figure  2.2  shows  a  branched  line 
system. 

Figure  2.2  shows  that  the  branched  line  model  asstimes 
that  the  pressure,  P,  at  the  junction  is  the  same  for  each 
line.  As  such,  the  load  impedance  for  the  input  line  (1) 
is  a  combination  of  the  end  impedances  for  lines  2  and  3. 
Thus,  Zli  is  now  given  by  (Kirshner  and  Katz,  1975:98)  : 

Z  =  .  (2-21) 

+  Zi3 
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III.  The  Compressor  Research  Facility  and  Model 

3.1  Introduction 

In  this  chapter,  the  geometry  of  the  CRF  and  the  scaled 
laboratory  model  are  detailed.  During  the  discussion, 
correlations  between  the  equations  presented  in  Chapter  II  and 
the  specific  geometry  of  the  facility  are  made. 

1^2  The  Compressor  Research  Facility 

As  illustrated  in  Figure  3.1,  the  CRF  is  69.8  ft  long  and 
includes  a  20  ft  diameter  inlet  section  (sf'ction  1),  a  10  ft 
diameter  flow  conditioning  section  (section  2),  and  a  3.3  ft 
diameter  section  (section  3)  at  the  inlet  of  the  test  article. 
Additionally,  a  2.5  ft  diameter  pipe  is  connected  to  the  inlet 
section  and  exhausted  outside  the  main  flow.  Mass  flow  rates 
through  the  facility  generally  range  from  20  Ibm/s  up  to  250 
Ibm/s . 

3.2.1  Section  Lengths.  The  lengths  of  sections  1  and  3 
(LI  and  L3)  can  vary  to  accommodate  different  con^ressors . 
Changes  in  LI  are  offset  by  an  equal  change  in  L3  such  that 
the  overall  facility  length  of  69.8  ft  remains  unchanged.  The 
maximum  and  minimum  lengths  of  both  sections  are  noted  in 
Figure  3.1  and  the  three  length  configurations  defined  for 
this  study  are  listed  in  Table  3.1. 


Table  3.1.  Length  Configuration  Definitions 


Configuration 

LI  (ft) 

L2  (ft) 

L3  (ft) 

Basic 

36.8 

19.8 

13.2 

Intermediate 

41.3 

19.8 

8.7 

Extended 

45.8 

19.8 

4.2 
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In  Table  3.1,  the  basic  length  configuration  is  defined 
such  that  LI  is  minimized  and  L3  is  maocimized.  The  extended 
length  configuration  is  defined  for  L3  maximized  and  LI 
minimized.  The  intermediate  length  configuration  is  defined 
for  both  LI  and  L3  between  their  maximum  and  minimum  lengths . 
For  all  configurations,  the  sum  of  LI,  L2,  and  L3  is  69.8  ft. 


3.2.2  Flow  conditioning  elements.  Section  2  contains  two 
sets  of  grids,  screens,  and  honeycomb,  called  flow 
conditioning  elements,  that  remove  axial  and  radial  turbulence 
in  the  throughflow.  These  elements  improve  the  quality  of  the 
throughflow  at  the  compressor  face  but  also  add  iitpedance  to 
the  system.  Figure  3.2  shows  a  breakdown  of  the  flow 
conditioning  elements  and  their  locations.  The  upstream  end, 
denoted  as  line  A-A  in  Figure  3.2,  is  39.4  ft  from  the  inlet 
of  the  CRF  (for  the  basic  length  configuration) . 

The  impedance  presented  by  the  screens,  honeycomb,  and 
grid  elements  is  modeled  analytically  as  area  reductions  at 
the  locations  where  they  are  mounted.  These  areas  are 
determined  by  the  porosity  of  the  individual  elements.  For 
example,  the  equivalent  open  area  of  the  18 -mesh  screen  is 
calculated  as  (Rae  and  Pope,  1984:76) : 


K 


open  area 
total  area 


(  1  -  wd^) 


where,  for  the  18 -mesh  screen: 


(3-1) 
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w  -  18  wires/ inch 


d  =  0.017  inch  (dia) 


To  find  the  diameter  representing  the  equivalent  open  area: 


=  0.482  = 


D  =  6.94  ft 


JC  ( 10  ft)2 


In  the  analytical  model,  the  18-mesh  screen  is  represented  by 
a  line  of  diameter  6.94  ft  and  length  0.017  in. 

The  radius  changes  between  the  individual  sections  and 
the  flow  straightening  elements  requires  that  the  facility  be 


A 


Line  A  -  A  is  located  39.4  ft  from  the  inlet  valves  (basic  length  configuration) 

Figure  3.2.  Flow  Conditioning  Elements 
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treated  as  a  system  of  cascaded  lines  in  the  fluid 
transmission  line  equations.  Thus,  the  gain  for  each  line  is 
calculated  and  multiplied  by  the  gains  of  the  other  lines  to 
find  the  overall  gain  for  the  facility. 

3.2.3 _ Pressure  Equalization _ EJJ2&4.  The  pressure 

equalization  pipe  is  detailed  in  Figure  3.3.  The  function  of 
the  pipe  is  to  provide  inlet  flow  to  a  large  cavity 
surrounding  the  outside  structure  of  sections  2  and  3  of  the 
facility.  Because  the  cavity  surrounds  the  outside  of 
sections  2  and  3,  it  is  not  in  the  main  flow  of  the  CRF.  The 


Figure  3.3.  Pressure  Equalization  Pipe 
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pipe  ensures  the  cavity  is  at  the  same  pressure  level  as 
section  1. 


For  the  purposes  of  this  study,  the  pressure  equalization 
pipe  is  simply  a  branched  line  that  draws  flow  from  the  20  ft 
diameter  section  and  exhausts  it  outside  the  main  flow.  From 
Figure  3.3,  it  is  clear  that  the  branched  line  is  conposed  of 
two  lines.  The  2.5  ft  diameter  main  line  is  conprised  of  the 
35.7  ft  and  17.8  ft  lines.  The  main  line  exhausts  to  the 
cavity  and  in  the  analytical  model  is  considered  a  branched 
line  (branching  from  section  1)  .  The  end  of  the  main  line  is 
treated  as  fully  open. 

The  48.0  ft  line,  also  of  2.5  ft  diameter,  branches  from 
the  main  line  to  a  valve.  In  the  analytical  model,  the  48.0 
ft  line  is  considered  a  branched  line  (branching  from  the  main 
line)  and  the  end  is  treated  as  fully  bloclced. 

3.2.4  Inlet  Valves.  The  five  inlet  butterfly  valves  are 
arranged  in  a  circular  pattern  about  the  centerline  of  the 
facility.  The  two  largest  valves  have  the  same  diameter,  36 
in,  while  the  other  three  valves  have  24  in,  20  in,  and  14  in 
diameters.  These  five  valves,  when  moved  while  the  facility 
is  operating,  produce  the  unsteady  disturbances  that  are  the 
input  signals  for  the  transmission  line  equations. 

For  the  analytical  model,  the  valves  are  not  considered 
individually.  Instead,  it  is  assumed  that  a  single  sinusoidal 


disturbance  emanates  from  the  point  vdiere  the  valves  are 
located  in  the  line  (inlet,  x  =  0) .  This  approximation  was 
made  in  order  to  sinplify  the  analysis;  because,  modeling  the 
individual  valves  would  be  an  effort  in  itself. 

The  transmission  line  equations  depend  not  only  on  the 
inlet  conditions  but  also  on  the  end  conditions.  Here,  too, 
the  geometry  is  approximated  to  simplify  the  analysis. 

3.2.5  Exit  Geometry.  The  CRF  exhausts  the  main  flow 
through  the  compressor  into  a  long,  taurus-shaped  pipe  and 
eventually  returns  the  flow  to  atmospheric  conditions.  For 
the  purposes  of  this  study;  however,  the  conpressor,  taurus- 
shaped  pipe,  and  eventual  return  of  the  flow  to  atmospheric 
conditions  are  all  ignored.  The  validity  of  this  assumption 
is  discussed  in  Chapters  VI,  VII,  and  VIII. 

Including  the  conpressor  in  the  analysis  requires 
detailed  knowledge  of  the  energy  added  to  the  throughflow  by 
each  test  article.  Since  each  test  article  behaves 
differently,  it  would  be  a  significant  task  to  find  a  model 
that  covers  all  of  the  potential  conpressors  to  be  used  in  the 
facility.  Because  of  this,  the  conpressor  and  all  downstream 
elements  are  ignored. 

The  analytical  model  instead  assumes  the  flow  in  section 
3  smoothly  necks  down  to  the  conpressor.  The  end  of  the 
transmission  line  is  considered  to  be  at  the  conpressor  face 
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and  the  compressor  face  is  assumed  to  be  at  the  end  of  section 
3,  69.8  ft  from  the  inlet  valves.  Modeling  the  system  exit  in 
this  manner  means  that  in  the  equation  for  end  iirqpedance  (2- 
15),  the  loss  coefficient,  Kj,,  has  a  constant  value  of  0.05 
(Kirshner  and  Katz,  1975:18).  While  far  from  exact,  this 
model  does  serve  the  purpose  at  hand,  which  is  to  find  the 
fundamental  frequencies . 


3.3  Theoretical  Limitations  on  the  CRF 

This  section  details  the  theoretical  limitations  on  the 
frequency  range.  Some  of  the  discussion  from  Chapter  II  that 
was  postponed  is  pic)ced  up  again  and  resolved. 

The  one-dimensional,  plane  wave  assumption  places  limits 
on  the  frequencies  for  which  the  Nichols  theory  is  valid. 
F.T.  Brown  (1962:548)  reports  that  for  a  valid  solution: 


f  «  ^  .  (3-2) 

Thus,  the  limiting  frequency  depends  on  the  largest  diameter 
line  in  the  system.  For  a  =  1130  ft/s  and  D  =  20  ft: 


f  « 


1130  ft/s 
WTt 


=  113  Hz  . 


(3-3) 


The  frequency  response  solution  is  limited  to  112  Hz  because 
of  this  restriction. 

Now  that  a  valid  frequency  range  has  been  established, 
recall  the  postponed  discussion  in  Section  2.1  on  the  number 
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of  resonant  peaks  for  single  constant  diameter  closed  and  open 
pipes  in  a  given  frequency  range.  The  statement  was  made  that 
under  certain  length  and  frequency  conditions,  the  blocked 
pipe  can  have  more  resonant  frequencies  than  the  open  pipe. 

At  this  time,  the  length  of  the  CRF  has  been  presented  as  68.9 

! 

I 

ft,  the  valid  frequency  range  as  112  Hz,  and  the  speed  of 
sound  can  be  assigned  the  standard  sea  level  value  of  1130.5 
ft/s.  Using  the  basic  pipe  resonance  equations  (2-1  and  2-2), 
the  number  of  peaks  for  siitple  blocked  and  open  lines  of 
length  69.8  ft  can  be  determined. 

For  the  open  line,  the  highest  resonant  frequency  below 
112  Hz  occurs  for  n  =  13  at  105.3  Hz.  This  means  an  open  line 
of  length  69.8  ft  will  have  13  resonant  frequencies  in  the 
valid  frequency  range.  On  the  other  hand,  for  the  blocked 
line  of  length  69.8  ft,  the  highest  resonant  frequency  below 
112  Hz  occurs  when  n  =  27  at  a  frequency  of  109.3  Hz.  This 
means  that  the  blocked  pipe  of  the  specified  length  will  have 
14  resonant  frequencies  (recall  only  odd  n  are  present  for  the 
blocked  pipe)  in  the  valid  frequency  range. 

This  is  not  a  rigid  analysis  of  the  system  but  rather  an 
attenpt  to  gather  additional  background  information  to  help 
with  the  analysis  of  the  results  of  the  stuc^*^. 
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The  model  is  constructed  to  incorporate  as  many  of  the 
transmission  line  characteristics  of  the  CRF  (full-scale 
facility)  as  possible.  The  length  changes  of  LI  and  L3  are 
incorporated  in  the  model,  along  with  the  pressure 
equalization  pipe,  the  inlet  valves,  and  the  flow  conditioning 
elements.  The  end  of  the  model  is  terminated  with  a  short 
line  and  exhausted  to  the  atmosphere.  The  schematic  is 
equivalent  to  Figure  3.1,  with  the  exception  of  the  radii  and 
end  condition. 

The  model  is  constructed  from  high  density  polyethylene 
plastic  tubing.  This  is  a  rigid  material  that  allows  the 
input  signal  to  traverse  the  line  without  being  absorbed  by 
the  structure. 

3.4.1  Model  Length  and  Radii.  Table  3.2  illustrates  the 
model  dimensions  in  comparison  to  the  full-scale  facility.  As 
shown  in  the  table,  the  model  length  is  equal  to  the  length  of 
the  full-scale  CRF.  Each  section  is  maintained  on  a  1:1  ratio 
with  its  counterpart  in  the  CRF.  The  model  was  constructed  to 
allow  changing  Ll  and  L3  from  the  basic  to  the  intermediate 
and  extended  length  configurations,  although  only  the  basic 
length  configuration  was  tested.  In  order  for  the  model  to 
fit  in  the  provided  lab  space,  the  tubing  was  coiled  and 
contained  in  a  44  in  square  box. 
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Unlike  the  length,  the  radii  are  scaled  significantly. 
Sections  1  and  2  and  the  pressure  equalization  pipe  are  scaled 
1:120;  whereas  section  3  and  the  terminating  line  are  scaled 
1:106.  The  different  scales  are  a  result  of  standard  sizes 
for  commercially  available  tubing. 

3.4.2  Flow  Conditioning  Elements  for  the  Model.  The 
model  is  constructed  such  that  one  of  three  flow  conditioning 
sections  can  be  used  in  section  2.  Two  of  the  sections  are 
used  to  determine  if  the  analytical  model  for  the  flow 
conditioning  elements  is  valid  while  the  third  is  used  to 
verify  the  response  of  the  full-scale  CRF  analyses.  The  three 
sections  are  illustrated  in  Figure  3.4. 

The  purpose  of  the  washer  configuration  (top,  Figure 
3.4)  is  to  physically  model  the  analytical  representation  of 

Table  3.2.  CRF  and  Model  Dimensions  (Basic  Length) 


Dimension 

Model 

CRF 

LI  (Basic) 

36.8  ft 

36.8  ft 

L2 

19.8  ft 

19.8  ft 

L3  (Basic) 

13.2  ft 

13.2  ft 

R1 

1.0  in 

10.0  ft 

R2 

0.5  in 

5.0  ft 

R3 

0.1875  in 

1.65  ft 
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the  flow  conditioning  elements  that  is  used  in  the  analysis. 
Since  the  analysis  represents  only  area  reductions  at  the 
locations  of  these  elements,  each  element  in  the  analysis  can 
be  physically  described  as  a  washer-  Thus,  the  first  of  the 
three  flow  conditioning  sections  was  constructed  by  installing 
washers,  of  inner  radii  equivalent  to  the  area  reductions 
calculated,  in  the  locations  where  screens  would  normally  be 
placed.  The  grid  and  honeycomb  have  very  large  porosity, 
therefore  large  equivalent  radii,  and  were  not  modeled  in  this 
flow  conditioning  section. 
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The  second  flow  conditioning  section  constructed  (middle. 
Figure  3.4)  is  closely  related  to  the  first  in  that  only  the 
screens  are  modeled.  Here,  the  washers  are  not  present; 
instead,  screens  identical  to  the  screens  in  the  CRF  are 
installed  in  the  locations  where  the  washers  are  installed  for 
the  first  flow  conditioning  section  model.  Data  taken  with 
this  section  can  be  compared  to  data  with  the  washer  section 
to  determine  if  the  area  reduction  technique  properly  models 
the  screens . 

The  third  flow  conditioning  section  (bottom.  Figure  3.4) 
contains  both  screens  and  honeycomb  in  the  locations  where 
they  are  installed  in  the  full-scale  CRF.  This  model 
represents  the  most  realistic  configuration  of  the  CRF  and  is 
used  primarily  to  verify  the  analytical  results  for  the 
response  of  the  full-scale  facility. 

3.4.3  Pressure  Ecfualization  Pipe.  The  pressure 
equalization  pipe  is  included  in  the  model  by  branching  a  0.25 
in  diameter  line  off  of  the  inlet  of  section  1.  Note  that  the 
diameter  is  scaled  1:120  from  the  full-scale  value  of  2.5  ft. 
The  length,  as  with  the  other  sections  of  the  model,  is  not 
scaled. 

The  inlet  to  the  pressure  equalization  pipe  is  located 
0.33  in  downstream  of  the  inlet  of  the  model  and  is  detachable 
so  the  effect  of  its  presence  on  the  model  response  can  be 
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assessed.  Except  the  difference  in  radii,  the  analytical 
representation  for  the  model  pressure  equalization  pipe  is 
identical  to  the  analytical  representation  of  the  full-scale 
pressure  equalization  pipe. 

3.4.4  Model  Inlet  and  Exit.  The  inlet  of  the  model  is 
designed  to  represent  the  five  butterfly  valves.  There  are 
five  input  ports  for  accepting  either  the  sinusoidal  signal  or 
throughf]  w.  The  different  valve  sizes  of  the  CRF  inlet  are 
not  represented.  Instead,  the  inner  diameter  of  all  ports  is 
0.25  in.  The  ports  are  installed  in  a  2  in  diameter  pipe  with 
the  relative  spacing  about  the  centerline  equivalent  to  the 
relative  spacing  of  the  valves  in  the  CRF  inlet. 

The  physical  exit  of  the  model  is  slightly  different  from 
the  analytical  representation  of  the  full-scale  CRF  exit. 
Where  the  CRF  exit  is  analytically  modeled  as  a  smooth 
contraction  to  the  compressor  face,  the  laboratory  model  is 
physically  terminated  with  a  sharp  contraction  to  the  end 
radius .  Flow  through  the  model  is  exhausted  to  the  atmosphere 
through  the  exit  line;  whereas,  the  analytical  representation 
of  the  full-scale  CRF  exit  does  not  assume  the  flow  returns  to 
atmospheric  conditions. 

The  1.0  in  long  end  termination  of  the  model  is  very 
short  in  length  when  coitpared  to  the  rest  of  the  model  and  is 
referred  to  as  an  orifice-terminated  line.  This  type  of 
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termination  requires  a  different  calculation  for  the  loss 
coefficient  than  is  required  for  the  full-scale  analysis. 
Whereas  for  the  full-scale  analysis,  was  0.05  for  all  end 
radius  conditions,  the  model  Ki,  depends  highly  on  the  radius 
of  the  teinnination  since  the  flow  returns  to  atmospheric 
conditions  at  the  exit.  For  the  model,  then,  the  following 
equation,  given  by  Kirshner  and  Katz  (1975:15-31),  applies  for 
the  loss  coefficient  in  Equation  (2-15) : 


K^  =  1.5 


R3 


(3-4) 


In  Equation  (3-4),  f  is  the  Darcy  friction  factor  defined  for 
laminar  pipe  flow  as  (Gerhart  and  Gross,  1985:422)  ; 
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and  for  turbulent  pipe  flow  as: 


0.316 


’  turbulent 
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However,  only  the  real  part  of  end  inpedance  is  modeled  using 
Equation  (2-15) .  Recall  from  Chapter  II  that  inductance. 
Equation  (2-16),  should  be  included  when  the  line  length  is 
short.  For  the  model,  the  terminating  line  is  short,  but 
finite,  and  the  inductance  is  included  as  the  imaginary  part 
in  the  end  impedance. 

Several  end  radii  are  studied  with  the  model.  Under  the 
open  condition,  the  end  radius  is  equal  to  R3  (0.1875  in). 
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For  the  blocked  condition,  the  end  radius  is  zero.  All  other 
end  terminations  are  between  the  blocked  and  open  conditions. 
These  other  radii  are  ViR3  (0.04689  in),  (0.09375  in)  and 
MR3  (0.141  in) . 
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This  chapter  is  divided  into  two  main  parts.  First,  the 
experimental  apparatus  and  procedure  are  described  for  the 
frequency  response  tests.  Second,  the  transient  response  test 
set-up  and  procedure  are  discussed.  Details  on  equipment 
specifications,  such  as  manufacturer,  model,  and  serial 
numbers,  are  available  in  Appendix  C. 

4.2  Frequency  Response  Set-Up  and  Procedures. 

The  end  goal  of  the  frequency  response  tests  is  to 
determine  the  gain  and  phase  shift  of  an  input  signal  at 
various  locations  along  the  line.  Figure  4.1  is  an  overall 
schematic  of  the  experimental  apparatus  used  to  accomplish 
these  goals. 

Figure  4.1  shows  that  pressurized  air  (house  air)  is  fed 
to  a  pneumatic  signal  generator  that  produces  a  sinusoidal 
pressure  wave.  The  sinusoidal  wave  is  fed  through  one  of  the 
five  input  ports  at  the  model  inlet  and  allowed  to  traverse 
the  length  of  the  model.  At  four  locations  along  the  model 
length,  the  input  signal  is  sensed  by  dynamic  pressure 
transducers  (PS,  PR2,  PR3,  and  PR4) .  The  signals  from  these 
transducers  are  airplified  and  fed  to  a  digital  data 
acquisition  system.  The  acquisition  system  is  connected  to 
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Figure  4.1.  Experimental  ^^paratus 


a  personal  cc»nputer  that  displays  the  signals  and  allo%ir8 
manipulation  of  the  digital  data.  The  manipulated  data  give 
arrplitude  and  phase  shift  information  for  the  pressure  data 
collected  at  each  transducer  location  along  the  line.  The 
data  are  recorded  and  the  gain  and  phase  shift,  relative  to 
the  inlet  transducer,  are  coitputed. 

Two  types  of  experiments  were  conducted  to  determine  the 
frequency  response  of  the  model.  The  first  type  is  termed  "no 
throughf low" ,  indicating  that  the  mass  flow  in  the  model  is 
provided  entirely  by  the  pneumatic  signal  generator.  In  the 
no  throughf  low  experiments,  the  signal  was  supplied  to  one 
inlet  port  and  the  rest  of  the  ports  were  blocked.  These 
experiments  were  conducted  to  gather  baseline  data  so  the 
effect  of  throughflow,  end  inpedance,  and  the  flow 
conditioning  elements  could  be  determined. 

The  second  type  of  experiment  is  termed  "with 
throughflow"  and  indicates  that  mass  flow  was  provided  by  both 
the  pneumatic  signal  generator  and  steady  state  throughflow. 
For  the  experiments  with  throughflow,  all  the  input  ports 
accepted  flow;  one  accepted  the  pneumatic  signal  and  the  other 
four  accept  throughflow  from  the  house  air  supply.  These 
experiments  were  conducted  to  find  the  frequency  response  of 
the  model  under  conditions  similar  to  the  flow  conditions  in 
the  CRF. 
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The 


house  air  is  maintained  at  100  psig  by  a  large  compressor. 
This  is  the  air  supplied  to  both  the  pneumatic  signal 
generator  and  the  model  inlet . 

The  steady  throughflow  and  alternating  flow  from  the 
pneumatic  signal  generator  parallel  each  other  prior  to 
entering  the  inlet  of  the  model.  As  such,  the  pneumatic 
signal  generator  has  its  own  mass  flow  meter  and  pressure 
regulating  valve.  The  throughflow  has  a  separate  mass  flow 
detector  and  pressure  regulating  valve  (not  shown  in  Figure 
4.1)  . 

The  mass  flow  meter  for  the  pneumatic  signal  generator 
flow  is  an  electronic  device  that  uses  temperature 
measurements  to  calculate  the  flow  rate  with  a  high  degree  of 
accuracy.  The  mass  flow  for  throughflow  is  determined  from 
transducer  readings  on  either  side  of  an  orifice  place  using 
the  following  equation  (Gerhart  and  Gross,  1985:486) : 

where  Co  is  the  discharge  coefficient  for  the  orifice  and  AP 
is  the  pressure  drop  across  the  orifice  plate. 

Two  regulating  valves,  one  upstream  of  the  pneumatic 
signal  generator  and  another  for  throughflow  control,  are  used 
to  control  the  input  pressure  levels  to  well  below  100  psig. 
The  valve  typically  produces  a  small  signal,  on  the  order  of 
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5  X  10'^  psig  (peak-to-peak)  in  section  1.  For  experiments 
with  throughflow,  the  alternating  signal  from  the  pneumatic 
signal  generator  is  introduced  along  with  throughflow.  The 
amplitude  of  the  alternating  pressure  signal  remains  on  the 
order  of  5  X  10*^  psig  (peak-to-peak)  but  the  steady  pressure 
is  increased  from  0.0  psig  (no  throughflow)  to  around  5  psig 
in  section  1. 

4.2.2  Pneumatic  Signal  Generator.  The  pneumatic  signal 
generator  uses  high  voltage  applied  to  a  piezoelectric  disk  to 
produce  a  waveform.  The  waveform  shape  and  frequency  are 
determined  by  input  from  an  electronic  function  generator. 
The  pneumatic  signal  generator  is  capable  of  producing 
sinusoidal,  triangular,  and  square  waves  at  frequencies  as 
high  as  1000  Hz.  For  this  study,  only  sinusoidal  input 
signals  below  113  Hz  were  produced. 

These  waveforms  are  fairly  noisy  below  about  80  Hz. 
Above  80  Hz,  the  valve  produces  much  cleaner  signals.  Figures 
4.2a  and  4.2b  illustrate  the  waveform  from  the  valve  at  6.3  Hz 
and  93  Hz. 

The  input  waveform  is  inportant  because  this  is  the 
signal  that  traverses  the  model  and  is  to  be  sensed  by  the 
pressure  transducers.  If  the  signal  has  poor  characteristics 
to  begin  with,  it  is  probable  that  the  same  signal  will  be 
difficult  to  detect  at  the  end  of  the  model. 
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4.2.3  Transducers .  Five  transducers  are  located  along 
the  length  of  the  model.  Figure  4.3  is  a  schematic  of  the 
locations  of  the  five  transducers. 

All  of  the  transducers  in  the  model  are  capable  of 
sensing  dynamic  signals.  Each  has  a  differential  pressure 
port  that  can  be  attached  to  a  stea(^  state  pressure  port  in 
the  line  and  remove  the  mean  pressure  signal.  PS,  PR2,  PR3, 
and  PR4  are  the  main  sensors,  set  up  with  differential 
pressure  to  sense  the  unsteady  flow  only.  The  mean  pressure 
transducer  is  not  used  as  a  dynamic  sensor.  Instead,  for  this 
transducer,  the  differential  pressure  port  is  open  to 
atmosphere  and  the  mean  gage  pressure  in  section  1  is  sensed. 
The  mean  gage  pressures  in  sections  2  and  3  are  calculated 
based  on  the  section  1  mean  gage  pressure  and  friction  losses 
through  the  model . 

The  resonant  frequency  of  the  cfynamic  transducers  used  is 
85  KHz  (before  installation) ,  well  above  the  highest  pressure 
input  frequency  of  112  Hz.  These  transducers  have  a  5  psig 
range  and  output  voltages,  after  amplification,  between  0  and 
6  volts. 

The  transducers  are  extremely  accurate  and  are  able  to 
detect  very  low  anplitude  signals.  After  anplif ication,  most 
signals  detected  during  model  testing  were  in  the  millivolt 
and  microvolt  range.  These  voltages  correspond  to  pressure 
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an^litudes  between  0.02  psig  for  millivolt  signals  and  5  X  10'® 
psig  for  microvolt  signals. 

Figures  4.4a  an^"^  4.4b  illustrate  the  shape  of  the 
pressure  signals  at  the  transducer  locations  PS,  PR2,  PR3,  and 
PR4  for  the  6.3  Hz  and  93  Hz  input  signals,  respectively.  The 
figures  show  that  the  input  signals  of  Figures  4.2a  and  4.2b 
are  transmitted  down  the  model  and  detected  by  all  of  the 
transducers.  In  Figure  4.4a,  for  the  6.3  Hz  signal,  the 
signal  amplitude  and  shape  are  not  as  clearly  defined  as  in 
Figure  4.4b.  This  is  a  result  of  the  input  signal  at  6.3  Hz 
(Figure  4.2a)  containing  more  noise  than  at  the  higher 
frequency  of  93  Hz  (Figure  4.2b). 

4.2.4  Data  Acquisition  System.  The  system  used  for  data 
acquisition  is  a  16  channel  unit  with  12-bit  resolution.  Each 
channel  is  referred  to  the  same  internal  cloclc  and  therefore, 
all  channels  acquire  data  simultaneously.  The  data 
acquisition  system  includes  a  personal  computer  with  software 
that  allows  manipulation  of  the  dynamic  signals. 

Pressure  levels  and  phase  shift  are  easily  determined 
using  the  spectral  analysis  software  program  (Nicolet, 
1991:19-11,19-12)  provided  by  the  data  acquisition  system 
manufacturer.  This  program  converts  time  domain  data  to  the 
frequency  domain. 


4-9 


4-] 


Figures  4.5a  and  4.5b  illustrate  the  anplitude  of  the  6.3 
Hz  and  93  Hz  signals  presented  in  Figure  4.4  after  spectral 
analysis  has  been  performed.  Due  to  the  input  signal  having 
more  noise  in  the  low  frequency  case  (Figure  4.5a),  the 
amplitude  of  the  frequency  domain  "spike"  is  lower  for  the  6.3 
Hz  signal. 

For  the  frequency  response  experiments,  the  saitpling 
frequency  was  kept  constant  at  1  KHz,  well  above  the  214  Hz 
aliasing  frequency.  Each  saitpling  required  approximately  3 
seconds  to  gather  4000  points  per  channel.  Primarily,  these 
values  were  determined  by  trial  and  error,  chosen  because  the 
data  acquisition  time  was  minimized  without  losing  clarity  of 
the  signal. 


4.2.5  Frequency  Data  Acquisition  Procedure.  Prior  to 
running  the  first  test  of  the  day,  the  equipment  was  turned  on 
and  allowed  to  warm  up  for  one  hour.  After  warm-up,  the 
individual  test  procedure  began  with  recording  the  ambient 
conditions,  configuring  the  data  acquisition  system,  and 
setting  up  the  model  geometry. 

Testing  began  by  applying  pressure  to  the  pneumatic 
signal  generator,  applying  steady  inlet  flow  for  test 
conditions  requiring  throughflow,  setting  a  frequency  on  the 
electronic  function  generator,  and  allowing  the  system  to 
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settle  for  a  few  seconds.  When  the  system  settled,  the  data 
acquisition  system  was  commanded  to  take  readings. 

Once  the  data  was  acquired,  mass  flows  and  mean  pressures 
were  recorded.  The  dynamic  signals  for  PS,  PR2,  PR3,  and  PR4 
were  transformed  to  the  frequency  domain  using  the  spectral 
analysis  program  and  the  anplitude  and  phase  shift  of  each 
transducer  signal  was  recorded.  A  different  frequency  was  set 
on  the  electronic  function  generator,  the  system  allowed  to 
settle,  and  the  process  repeated  by  commanding  the  acquisition 
system  to  acquire  data. 

4.3  Transient  Experimental  Set -Up  and  Procedures. 

The  transient  tests  were  performed  using  a  solenoid  valve 
in  place  of  the  pneumatic  signal  generator.  The  procedures 
were  the  same  as  with  the  frequency  experiments  with  the 
exception  of  the  actual  data  acquisition. 

In  the  transient  tests  with  positive  step  inputs,  the 
system  was  maintained  at  atmospheric  pressure  until  the  valve 
was  opened.  The  acqui'^ition  system  was  triggered  just  prior 
to  opening  the  valve. 

Once  the  valve  opened,  the  acquisition  system  accjuired 
data  for  approximately  15  seconds.  After  the  data  acquisition 
was  complete,  the  valve  was  closed.  For  these  tests  the 
sampling  frequency  was  4.35  KHz  with  65,000  data  points 
acquired  per  channel . 
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For  negative  step  experiments,  the  solenoid  valve  was 
opened  and  the  model  allowed  approximately  one  minute  to  reach 
steady  state  pressure.  Just  prior  to  closing  the  valve,  the 
data  acquisition  system  was  triggered.  The  valve  was  closed 
and  data  acquired  for  15  seconds  at  4.35  KHz. 


V.  Cofflputer  Program 


5.1  Introduction 

The  program  calculates  gain  and  phase  at  discrete 
frequencies  by  stepping  through  the  valid  frequency  range  in 
large  increments.  After  the  entire  range  is  swept,  the 
resonant  frequencies  are  determined.  The  program  then 
restarts  the  sweep,  iterating  around  these  frequencies  in 
small  increments. 

5.2  Inputs  .and  Output 

The  program  accepts  several  inputs  and  produces  multiple 
output  variables.  The  inputs  include  geometry  variables, 
ambient  conditions,  and  mass  flow  rate. 

The  output  is  contained  in  several  arrays.  Each  array 
contains  either  the  gain  or  phase  shift  at  a  specified 
location  in  the  line  for  the  discrete  frequencies  in  which  the 
equations  are  solved. 

Ambient  conditions  are  used,  along  with  mass  flow  rate 
and  geometry,  to  solve  for  velocities  and  pressure  drops  in 
the  line.  The  velocities  in  each  section  are  calculated 
through : 


5-1 


and  the  static  pressure  through  Bernoulli's  equation  for 
incoirqpressible  flow: 

^  -  AP  .  (5-2) 

In  Equation  (5-2),  AP  is  (Gerhart  and  Gross,  1985:415): 

AP  -  £f{f  (5-3) 

for  a  smooth-walled  pipe  of  diameter,  D,  and  length,  L,  where 
f  is  defined  in  Equation  (3-5)  for  laminar  flow  and  Equation 
(3-6)  for  turbulent  flow.  Through  the  flow  conditioning 
elements,  AP  is  calculated  using  equations  developed  by 
Luchuclc  and  Priclcett  (1984:10,20): 

AP  =  (2.541  -  2.837P)  (5-4) 

where,  for  a  screen  and  honeycomb  set,  P  is  the  product  of  the 
porosities  of  the  individual  elements.  For  the  CRF  flow 
conditioning  elements,  P  =  0.06. 

The  speed  of  sound  is  calculated  as: 

a  =[  Y  r]Vi  (5-5) 


and  the  Mach  number  follows  as: 
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a 

The  Mach  number  is  used  in  the  equations  for  characteristic 
iitpedance  (2-17),  propagation  operator  (2-18),  and  line  gain 
(2-19)  . 

^>3 .Program  LimitatJLons 

The  program  is  restricted  both  by  theory  and  by  the 
programming  language.  The  range  of  frequencies  calculated  is 
limited  to  below  113  Hz,  due  to  the  restrictions  on  theory 
presented  in  Equation  (3-2). 

For  small  argxjments,  the  programming  language,  MATLAB", 
will  compute  the  Bessel  functions  in  the  equations  for 
impedance  (2-7)  and  admittance  (2-8) .  These  intrinsic 
functions  wor)c  well  for  the  solution  of  the  scaled  laboratory 
model;  however,  due  to  the  round-off  errors  associated  with 
the  CRF  dimensions,  these  functions  do  not  yield  an  answer  for 
the  full-scale  solution.  When  full-scale  dimensions  are 
entered,  the  program  uses  an  algebraic  approximation  devised 
tjy  F.T.  Brown  (Kirshner  and  Katz,  1975:83)  to  solve  for  Z  and 
Y.  Although  the  intrinsic  functions  are  used  when  possible, 
the  approximation  is  generally  twice  as  fast  and  nearly  as 
accurate  as  the  intrinsic  functions. 
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A  wealth  of  data  already  exist  on  fluid  transmission 
lines  as  a  result  of  previous  work.  These  data  are  well 
established  and  were  used  to  verify  the  program.  Table  5.1 
lists  the  conditions  of  verification  and  sources  used  to 
validate  the  output  from  this  conputer  program.  All  of  the 
configurations  were  tested  using  both  the  intrinsic  Bessel 
functions  and  Brown's  approximation.  Although  the  data  is  not 
presented  here,  the  results  of  Runs  1.1  through  1.6  in  Table 
5.1  showed  that  this  computer  program  yields  the  same  solution 
as  presented  in  the  sources,  regardless  of  the  calculation 
method  used  for  the  Bessel  functions. 

Additional  program  verification  was  acconplished  using 
the  laboratory  model  of  the  CRF.  In  Chapter  VII,  the  results 
from  the  model  tests  are  reviewed  in  light  of  the  analytical 
results  obtained  for  the  full-scale  CRF.  First,  the  full- 
scale  results  obtained  from  several  runs  of  the  conputer 
program  are  presented  in  the  next  chapter. 
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Table  5.1.  Verification  Matrix  for  Conqputer  Program 


Coit4)uter 

Run  Nxamber 

Configuration 

Tested 

Con^ared  against 
published  data 

1.1 

Attenuation 

Karam  and  Franke, 
1967:375 

1.2 

Volume  Terminated 

Franke  et  al . , 
1969:413 

1.3 

Cascaded 

Franke,  Karam, 
and  Lymburner, 
1970:9 

1.4 

Branched 

Franke , 

Malanowski ,  and 
Martin,  1972:18 

1.5 

Orifice- terminated 

Franke,  Karam, 
and  Lymburner, 
1970:12 

1.6 

Mean  flow 

Franke,  Karam, 
and  Lymburner, 
1970:12 
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I  Scale  CRF  Results 


6.1  Introduction 

In  this  chapter,  the  results  from  several  ccmsputer  runs 
are  presented  for  the  full-scale  facility.  Several  flow  and 
geometry  configurations  are  investigated  and  the  relative 
importance  of  each  are  discussed. 

All  of  the  results  presented  in  this  chapter  are  for  the 
full-scale  facility  with  the  flow  straighteners  and  pressure 
equalization  pipe.  The  lengths  of  sections  1  and  3  are 
varied,  as  are  the  mass  flow  rates  and  compressor  radius. 

£.2  Tabulated  Results 

Tables  6.1,  6.2,  and  6.3  give  the  results  of  63  computer 
runs  for  varied  configurations  of  the  facility.  Each  table 
represents  specific  length  geometry.  The  basic  length 
configuration  is  represented  in  Table  6.1,  the  intermediate 
length  configuration  in  Table  6.2,  and  the  extended  length 
configuration  in  Table  6.3.  For  all  runs,  the  ambient 
pressure  is  14.3  psig  and  the  ambient  temperature  is  72  **F, 
resulting  in  the  same  speed  of  sound  for  all  cases. 

In  the  tables,  the  end  impedance  is  calculated  through 
Equation  (2-15)  except  for  the  no  throughflow  runs  (8.18  - 
8.21,  9.18  -  9.21,  and  10.18  -  10.21).  In  the  no  throughflow 
runs,  the  end  termination  is  artificially  assigned  a  unit 
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Table  6.1.  CRF  Basic  Length  Results** 
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Table  6.2.  CRF  Intermediate  Length  Results* 
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Table  6.3.  CRF  Extended  Length  Results 
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length  of  1.0  ft  and  the  end  impedance  is  calculated  as  the 
input  impedance  of  the  1.0  ft  line  using  Equation  (2-20). 
These  runs  required  this  end  impedance  calculation  because  the 
zero  mass  flow  rate  condition  forces  Z^,  -  0  in  Equation  (2- 
15) ,  causing  the  solution  to  Ise  undefined  when  Equation  (2-15) 
is  used. 

6.2.1  Fundamental  Frequency.  The  most  notable  difference 
between  the  three  tables  is  the  different  fundamental 
frequencies  for  each  length  configuration.  Specifically/  note 
that  the  frequency  varies  frc»B  a  low  of  5.5  Hz  (run  10.1/ 
Table  6.3)  to  a  high  of  6.9  Hz  (run  8.6  -  8.17/  8.20/  8.21/ 
Table  6.1)  depending  on  the  length  configuration  and  the  end 
impedance . 

6.2.2  Length  Configuration  Dependence.  For  all  mass  flow 
and  end  impedance  conditions/  the  extended  configuration 
yields  the  lowest  fundamental  frequencies.  For  the  blocked 
cases  (runs  8.1/  9.1/  and  10.1)/  the  frequency  drops  by  0.5  Hz 
from  the  basic  to  the  intermediate  configuration  and  by 
another  0.5  Hz  from  the  intermediate  to  the  extended 
configuration.  The  lowest  drop  in  fundamental  frequency  is 
noted  at  the  more  open  conditions  where  the  compressor  radius 
is  equal  to  or  larger  than  ^3.  Under  these  conditions/  the 
frequency  drops  0.5  Hz  from  the  basic  to  intermediate  length 
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configuration  and  0.3  Hz  fr<m  the  interziediate  to  extended 
configuration . 

In  all  cases /  the  fundamental  frequency  is  lowest  when  LI 
is  maximized  at  45.8  ft.  As  a  check  on  these  results, 
consider  the  basic  pipe  resonance  equations  [(2-1)  and  (2-2)] . 
Recall  that  for  the  simple  lines  represented  in  Equations  (2- 
1)  and  (2-2),  longer  lines  have  lower  fundamental  frequencies, 
provided  the  speed  of  sound  is  constant.  This  is  in 
accordance  to  what  is  found  by  comparing  data  in  the  tables. 
Although  the  overall  length  is  constant  at  69.8  ft,  LI  drives 
the  frequency  low  when  it  is  maximized,  similar  to  the  basic 
pipe  where  the  long  pipe  drives  the  fundamental  frequency  low. 

The  reason  for  this  dependence  on  LI  and  L3  is  that  the 
input  signal  partially  reflects  off  of  the  junction  between 
sections  1  and  2,  where  an  impedance  mismatch  occurs.  The 
impedance  mismatch,  caused  by  the  different  dicuneters  of 
sections  1  and  2,  is  located  a  different  distance  from  the 
inlet  for  the  three  cases.  Signals  have  farther  to  travel  in 
the  intermediate  and  extended  configurations  before  partial 
reflection,  resulting  in  lower  fundamental  frequencies. 

6.2.3  End  Impedance  Dependence.  In  general,  higher  end 
impedance  (a  small  compressor  and/or  a  high  flow  rate)  results 
in  lower  fundamental  frequencies.  This  is  expected  from  the 
basic  pipe  resonance  equations  [(2-1)  and  (2-2)],  where  the 
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blocked  (infinitely  large  end  in^dance)  fundamental  frequency 
is  one  half  of  the  open  (zero  end  impedance)  fundamental 
frequency. 

From  Tables  6.1,  6.2,  and  6.3,  it  is  clear  that  the  end 
impedance  is  higher  for  small  compressors  and  high  mass  flow 
rates,  as  would  be  expected  from  Equation  (2-15) .  In  fact, 
for  the  MR3  cases  (Runs  8. 2-8. 5,  9. 2-9. 5,  10.2-10.5),  the 
combination  of  a  small  compressor  and  high  flow  rates  causes 
the  fundamental  frequency  to  approach  that  of  the  blocked  line 
(Runs  8.1,  9.1,  and  10.1).  The  result  of  the  high  end 
impedance  is  the  fundamental  frequency  is  low,  close  to  the 
frequencies  found  for  the  blocked  line  cases. 

Comparing  the  unblocked,  no  throughflow  cases  (8.18  - 
8.21,  9.18  -  9.21,  10.18  -  10.21)  with  the  blocked  line 
response  (8.1,  9.1,  10.1),  it  is  evident  that  increasing  the 
end  radius  from  zero  changes  the  fundamental  frequency.  This 
is  due  to  the  fact  that  for  the  blocked  line  the  portion  of 
the  input  signal  that  propagates  to  the  end  is  reflected 
entirely.  For  partially  or  fully  open  ends;  however,  only  a 
portion  of  the  signal  is  reflected  at  the  exit.  Part  of  the 
signal  propagates  through  the  opening  and  exits  the  line. 
This  effectively  changes  the  transmission  line  characteristics 
and  the  result,  for  the  no  throughflow  conditions,  is  a  change 
in  fundamental  frequency  from  the  blocked  line  case. 
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with  throughflow  in  the  line,  wave  propagation  is  again 
affected,  now  due  to  the  interaction  of  the  throughflow 
velocities  with  the  propagating  signal.  These  interactions 
cause  changes  in  the  characteristics  of  the  line  without 
throughflow  (Katz,  Housner,  and  Eisenberg:  1972,  272), 
resulting  in  different  resonant  frequencies  for  the  cases  with 
throughflow.  Runs  with  throughflow  (8.2  -  8.17,  9.2  -  9.17, 
10.2  -  10.17)  can  be  coit^ared  to  runs  without  throughflow 
(8.18  -  8.21,  9.18  -  9.21,  10.18  -  10.21)  to  illustrate  this 
effect.  For  the  higher  mass  flow  rate  runs  where  ih  ^  180 
Ibm/s,  the  throughflow  velocity  in  section  3  can  be 
significant  coir^jared  to  the  speed  of  sound.  At  high  flow 
rates,  velocities  at  the  compressor  face  can  be  as  high  as  500 
ft/s  (M  «  0.5)  .  These  high  velocities  affect  the  propagation 
of  the  acoustic  signal  by  altering  the  velocity  of  the  wave. 
Waves  propagating  in  the  direction  of  throughflow  travel  at 
the  speed  of  sound  plus  throughflow  velocity;  whereas,  waves 
propagating  against  the  throughflow  direction  have  a  velocity 
of  the  speed  of  sound  minus  the  throughflow  velocity.  Because 
of  the  different  velocities,  the  interaction  between  the  waves 
propagating  in  the  two  directions  is  different  than  when  no 
throughflow  conditions  exist. 

6.2.4  Fundamental  Pealcs .  The  anplitude  of  the 
fundamental  pea)cs  is  highest  for  high  end  iitpedance  because 
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most  of  the  waves  are  reflected  at  the  end  and  not  allowed  to 
propagate  out  of  the  line.  For  example/  note  that  for  all  of 
the  blocked  conditions  (Runs  8.1,  9.1,  and  10.1)  the 
fundamental  peak  is  quite  large,  around  70  dB.  These  high 
peaks  are  a  result  of  the  large  radius  and  blocked  end. 
Additional  details  on  peak  levels  are  available  in  J^pendix 
B. 

For  the  fully  open  cases  with  no  mass  flow  (8.21,  9.21, 
and  10.21),  the  fundamental  peak  drops  to  just  above  half  the 
decibel  level  of  the  blocked  cases.  The  primary  reason  for 
this  drop  in  amplitude  is  that  the  open  condition  provides 
very  low  impedance.  Much  of  the  signal  is  propagated  out  of 
the  line,  leaving  little  to  be  reflected. 

In  cases  with  throughflow,  the  amplitudes  of  the 
fundamental  peaks  are  reduced  not  just  below  the  blocked 
condition,  but  also  below  the  open  condition  (Runs  8.14  - 
8.17,  9.14  -  9.17,  and  10.14  -  10.17).  The  peaks  are  highest 
for  high  mass  flows  and  small  compressor  radii .  This 
indicates  that  throughflow  velocities,  when  added  to  the  speed 
of  the  propagating  waves,  aids  the  signals  in  propagating  out 
of  the  system;  thereby  reducing  the  gains. 

6.3  Graphical  Results 

Six  figures  are  presented  and  analyzed  in  this  section 
for  the  valid  frequency  range.  The  solution  for  the  entire 
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frequency  spectrum  is  presented  for  the  blocked  line  response 
of  the  basic,  intermediate,  and  extended  length 
configurations.  Further,  three  plots,  with  the  con?)ressor 
radius  fixed  at  are  presented  for  a  single  mass  flow  rate 
and  varied  length  configurations. 

The  graphs  present  the  overall  gain  and  the  phase  shift 
of  the  signal  at  the  end  of  the  CRF  as  a  function  of 
frequency.  The  phase  shift  was  forced  to  limits  of  ±  180*  by 
using  an  intrinsic  function  in  MATLAB™. 

6.3.1  Blocked  Line  Results.  Figures  6.1  through  6.3 
illustrate  the  blocked  line  response  for  the  basic, 
intermediate,  and  extended  length  configurations, 
respectively.  These  curves  represent  hypothetical  situations 
because,  with  a  conpressor  installed  and  under  testing,  the 
facility  is  definitely  open  to  some  degree.  Although  they 
represent  a  hypothetical  situation,  these  conditions  are  a 
worst  case  scenario  since  the  fundamental  frequencies  are  the 
lowest  of  any  case  studied  (Tables  6.1  through  6.3)  and  they 
have  the  highest  fundamental  peak  airplitudes.  Both  low 
fundamental  frequencies  and  high  fundamental  peak  anplitudes 
present  problems  for  the  control  system  design  (Hull; 1993) . 

The  peaks  represent  the  gain  of  the  CRF  at  the  resonant 
frequencies .  The  lowest  frequency  peak  is  the  fundamental  and 
the  rest  are  harmonics.  There  are  14  peaks  present  in  all 
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Run  8.1 


Run  9.1 


plots,  indicating  the  blocked  behavior  expected  (Chapter  III) 
fr<xn  the  basic  pipe  resonance  equations  (2-1)  . 

Note  that  the  higher  frequency  peaks  are  in  the  same 
decibel  range  as  the  fundamental.  All  of  the  peaks  are  high 
and  very  sharp.  In  contrast,  off-peak  gains  are  low  and  more 
rounded.  The  corresponding  phase  plots  show  the  phase  angle 
crossing  zero  at  the  peak  frequencies  with  a  nearly  infinite 
slope,  further  illustrating  the  sharpness  of  the  peaks.  The 
off-peak  phase  angles  are  a  maximum  at  ±180  degrees  as  is 
expected  for  low  gains. 

Note  also  the  curves  are  generally  the  same  shape  for  all 
three  length  configurations.  However,  note  how  drastically 
different  the  resonant  frequencies  are  for  the  different 
length  configurations  above  20  Hz.  Above  20  Hz,  the  resonant 
peaks  clearly  occur  at  different  frequencies  for  the  three 
length  configurations,  indicating  the  entire  frequency 
spectrum  is  affected  by  LI  and  L3. 

6.3.2  Results  for  a  Fixed  Radius _ Compressor.  The 

remaining  three  figures  show  the  response  for  an  exit  radius 
0.825  ft  (^3)  and  a  mass  flow  rate  of  180  Ibm/s.  Figures 
6.4,  6.5,  and  6.6  show  the  response  for  the  basic, 
intermediate,  and  extended  length  configurations, 
respectively.  The  curves  for  a  compressor  radius  of  ^4R3  and 
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Figure  6.5.  CRF  Frequency  Response  for  W3  End  Radius 
Conditions:  R^  =  0.825  ft,  ih  =  180  Ibm/s.  Fundamental  Peak: 


Run  10.8 


mass  flow  rates  of  20  Ibm/s,  100  Ibrn/s,  and  250  Ibm/s  are 
located  in  Appendix  A,  Figures  A.l  -  A. 9. 


Comparing  Figures  6.4,  6.5,  and  6.6,  it  is  clear  that  the 
resonant  frequencies  are  dependent  on  the  length 
configuration.  This  is  the  same  result  as  with  the  blocked 
line  plots  in  Figures  6.1  -  6.3.  For  all  cases  in  Figures  6.4 
-  6.6,  the  number  of  peaks  in  the  spectrum  is  lower  than  for 
the  blocked  case,  varying  only  slightly  between  12  and  13. 
The  lower  number  of  peaks  is  expected  since  the  line  is  no 
longer  blocked  (Chapter  III  discussion) . 

All  figures  exhibit  the  Scune  peak  levels,  relative  to  the 
fundamental,  as  in  the  blocked  cases.  The  gains  are  lower 
than  for  the  blocked  cases,  but  rise  with  increasing  mass  flow 
conditions  (Appendix  A)  .  The  lower  gains  are  again  a  result  of 
the  wave  propagating  out  of  the  system  through  the  open  end. 
This  also  results  in  smoother  phase  transitions  through  zero 
and  increased  spacing  between  the  peaks,  as  compared  to  the 
blocked  line  curves . 


In  summary,  the  data  presented  yields  a  fundamental 
frequency  for  the  CRF  between  5.5  Hz  and  6.9  Hz,  depending  on 
the  length  configuration  and  end  impedance  conditions.  The 
gains  for  the  CRF  range  from  70  dB,  for  the  blocked  cases, 
down  to  13  dB  for  the  open  end  conditions  with  mass  flow.  The 
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frequency  trends  and  number  of  peaks  in  the  spectrum  checks 
with  the  trends  expected  from  the  basic  pipe  resonance 
equations  at  varied  length  and  end  impedance  conditions. 

However,  the  gains  cannot  be  checked  against  any  basic 
equations  since  none  exist.  It  turns  out  that  when  the  gains 
are  used  to  determine  the  damping  characteristics  of  the  CRF, 
the  amplitudes  indicate  a  very  lightly  damped  system 
(Hull: 1993) .  This  presents  a  problem  for  the  control  system 
designer.  Further  discussion  on  this  problem  is  postponed 
until  after  the  experimental  results  are  presented. 
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VII.  Experimental  Rftaulta 


T.LJatrQductiQn 

In  this  chapter,  the  results  of  several  experiments  with 
the  CRF  model  are  presented  and  discussed.  The  first 
discussion  centers  on  gain  and  phase  shift  results  from  data 
taken  at  discrete  frequencies  for  no  throughflow  conditions. 
These  results  are  compared  with  analytical  solutions  for  the 
model  at  the  flow  and  geometry  conditions  of  the  individual 
experiments .  The  second  section  discusses  results  for 
experiments  with  throughflow.  In  the  last  section,  transient 
response  test  results  for  the  model  are  presented  and 
analyzed.  All  results  presented  are  for  the  basic  length 
configuration. 

7.2  No  Throughflow  Model  Response 

Several  data  sets  without  throughflow  were  taken  using 
the  pneumatic  signal  generator  as  a  sinusoidal  pressure  input . 
Only  a  sampling  of  the  results  is  presented  here;  additional 
plots  are  located  in  Appendix  A. 

The  condition  of  no  throughflow  simply  means  that  the 
mass  flow  in  the  model  is  provided  only  from  the  sinusoidal 
input.  As  such,  mean  mass  flow  rates  for  these  experiments 
are  extremely  low,  on  the  order  of  2  X  10’*  slug/s.  This,  in 
turn,  yields  velocities  on  the  order  of  1  X  10'*  ft/s  in 
section  3  of  the  model.  Arguably,  these  conditions  are  not 
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representative  of  the  realistic  flow  in  the  CRF  where  the  mass 
flow  rate  and  velocities  are  much  higher.  However,  these 
results  provide  insight  into  the  accuracy  of  the  analysis  of 
the  full-scale  facility  and  are  presented  for  that  purpose. 

Two  groups  of  experiments  are  discussed  in  the  following 
sections.  The  first  group  is  comprised  of  four  runs  for  the 
model  configured  with  screens  and  honeycomb  in  the  flow 
conditioning  section.  The  end  radius,  Re,  varies  with  each 
run.  These  results  are  compared  with  full-scale  analyses 
for  the  same  flow,  length,  and  end  conditions. 

The  second  group  of  experiments  is  nearly  identical  to 
the  first,  except  for  the  flow  conditioning  elements  in  the 
model.  The  second  set  of  data  are  taken  with  the  washers 
installed  in  section  2  instead  of  the  screens  and  honeycomb. 
These  results  are  compared  to  the  model  results  in  the  first 
group  to  determine  the  effects  of  modeling  the  flow 
conditioning  elemeni.s  as  area  reductions. 

7.2.1  Screens  and  Honeycomb:  Blocked  Line  Response.  The 
blocked  line  response  for  the  basic  length  model  with  screens 
and  honeycomb  in  the  flow  conditioning  section  is  presented  in 
Figures  7 . la  through  7 . Ic^ .  Data  points  are  represented  by 

^Figure  numbers  appended  with  "a"  show  results  for  the 
end  transducer  (PR4)  as  compared  to  the  inlet  transducer  (PS)  . 
Figures  appended  with  "b"  and  "c"  show  results  for  the 
transducers  near  flow  conditioning  elements  (PR3  and  PR2, 
respectively)  .  Transducer  locations  are  defined  in  Figure  4.3. 
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igure  7.1a.  Model  Blocked  Line  Frequency  Response  (PR4) .  Conditions 
slug/s,  basic  length,  screens  and  honeycomb.  Fundamental  peak;  6.3 
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Model  Blocked  Line  Frequency  Response  (PR3).  Conditions:  R, 
.c  length,  screens  and  honeycomb.  Fundamental  peak:  6.3  Hz. 
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Model  Blocked  Line  Frequency  Response  (PR2) .  Conditions 
lasic  length,  screens  and  honeycoirib.  Fundamental  peak:  6.3 


circles  and  the  analytical  solution  for  the  model  is 
represented  by  the  solid  line.  The  full-scale  results  for  the 
same  flow  conditions  are  illustrated  in  Figures  7.2a  through 
7.2c.  There  is  no  full-scale  experimental  data;  therefore, 
only  the  analytical  curves  are  presented  in  full-scale  plots. 

The  fundamental  frequency  for  the  model  is  6.3  Hz; 
whereas,  for  the  full-scale  solution,  the  fundamental 
frequency  occurs  at  6.5  Hz.  The  fundamental  frequencies  of 
the  model  and  CRF  are  nearly  the  same,  as  would  be  expected 
from  the  radius  scale  study  performed  (Appendix  B) . 

Reviewing  the  gains  in  Figures  7.1a  -  7.1c,  note  that 
most  of  the  data  fall  on  the  analytical  curves  quite  well  at 
frequencies  below  50  Hz.  Above  50  Hz,  the  data  fall  in  the 
vicinity  of  the  curves  until  around  90  Hz  where  the  data  and 
analyses  diverge.  The  corresponding  phase  plots  show  the  same 
behavior  although,  in  general,  the  phase  data  follow  the 
analyses  more  closely  than  the  gain  data,  especially  at  higher 
frequencies.  The  high  frequency  discrepancy  between  theory 
and  experiment  is  common  throughout  all  of  the  plots. 
Although  the  cause  of  the  discrepancy  is  unknown,  it  is  most 
likely  due  to  losses  in  the  model  that  are  not  accounted  for 
in  the  analysis,  such  as  the  curvature  of  the  coiled  tubing. 

The  gains  vary  significantly  for  the  different  transducer 
locations.  While  PR4  (Figure  7.1a)  decibel  levels  are  all 
greater  than  zero,  both  PR3  (Figure  7 . lb)  and  PR2  (Figure 
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CRF  Blocked  Line  Frequency  Response  (PR4) .  Conditions 
length.  Fundamental  peak:  6.5  Hz. 
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7.1c)  gains  fall  below  0  dB  at  some  frequencies.  This  is 
simply  an  indication  that  gains  are  a  function  of  location  in 
rhe  line.  The  basic  behavior  of  the  line  is  not  dependent  on 
location  since  the  peaks  and  valleys  occur  at  the  same 
frequencies  for  all  three  locations. 

Finally,  comparing  Figures  7.1a  -  7.1c  with  7.2a  -  7.2c, 
it  is  obvious  that  the  gain  decibel  levels  are  significantly 
higher  for  the  full-scale  analysis  than  for  the  model  analysis 
and  data.  This  is  a  result  of  radius  scaling.  Additional 
details  on  radius  scaling  from  a  20  ft  diameter  line  to  a  2  in 
diameter  line  and  the  associated  peak  decibel  levels  are 
located  in  Appendix  B. 

Four  important  characteristics  of  the  model  are  found 
from  the  blocked  line  analysis,  where  infinite  exit  impedance 
exists.  First,  the  data  and  theory  are  in  good  agreement  for 
this  configuration  below  90  Hz.  Second,  the  amplitude  of  the 
resonant  peaks  are  dependent  on  location  in  the  line.  Third, 
the  model  peaks  are  significantly  lower  than  the  full-scale 
peaks  for  the  same  flow  and  geometry  configurations.  Fourth, 
both  the  theoretical  and  experimental  results  place  the 
fundamental  frequency  for  the  model  at  6.3  Hz.  This  is  very 
close  to  the  full-scale  fundamental  frequency  of  6.5  Hz, 
illustrating  that  the  experiment  verifies  the  full-scale, 
blocked  line  frequency  results. 
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7.2.2  Screens  and  Honeycoiab^  R-  "  Figures  7.3a 
through  7.3c  are  presented  in  this  section  for  the  basic 
length  configuration  with  screens  and  honeycomb  as  flow 
conditioning  elements  and  the  end  radius  fixed  at  ^3.  The 
mass  flow  is  again  due  only  to  the  sinusoidal  input  from  the 
pneumatic  signal  generator  and  for  these  figures,  is  computed 
to  be  2.2  X  10'*  slug/s. 

In  Figures  7.3a  through  7.3c,  again  note  the  data  points 
fall  closely  to  the  analytical  curves.  However,  in  this  case, 
the  agreement  is  not  as  good  as  in  the  blocked  case.  This  is 
a  result  of  the  end  impedance  equation  used  for  the  analysis. 
For  blocked  lines,  the  end  impedance  is  infinite  and  no 
approximation  of  its  value  is  necessary.  However,  for  an 
orifice-terminated  line,  the  end  impedance  is  not  well 
understood  and  approximations  for  the  value  must  be  made 
[Equation  (3-4)].  The  end  impedance  equation  used  (2-15)  for 
the  model  analyses  was  determined  by  trial  and  error.  Other 
representations  in  the  analysis  were  tried  but  were  much  less 
successful.  One  representation  tried  was  to  use  the  Equation 
(4-1)  for  the  pressure  drop  across  and  orifice  plate.  This 
representation  was  used  by  Ross  (1983:20)  with  successful 
results .  A  second  representation  tried  was  to  treat  the 
orifice  as  a  short,  open  line. 

Computing  the  frequency  response  using  Equation  (2-15) ; 
however,  produced  the  best  match  between  analytical  and 
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Model  No  Throughflow  Frequency  Response  {PR3)  for  R^  =  HR3 .  Conditions: 
rh  =  2.2  X  10'®  slug/s,  basic  length,  screens  and  honeycomb. 


theoretical  results.  The  use  of  different  end  Impedance 
equations,  in  general,  vertically  shifted  the  theoretical 
curve  and  slightly  changed  its  shape.  It  is  possible  to 
devise  a  better  analytical  representation  of  the  end  impedance 
by  further  trial  and  error. 

Comparing  again  the  blocked  and  ^3  cases,  note  that  the 
peak  gains  are  reduced  as  the  end  condition  changes  from 
blocked  to  partially  open.  Particularly,  the  fundamental  peak 
for  PR4/PS  (Figure  7.3a)  drops  to  20  dB  from  the  blocked  line 
level  of  30  dB  (Figure  7.1a).  From  the  analysis  in  Chapter 
VI,  this  is  an  expected  response  when  the  end  radius  is 
increased.  The  peaks  for  the  full-scale  CRF  also  drop  in 
magnitude  with  the  increased  end  radius.  These  curves  are  in 
Appendix  A  (Figures  A.  10  -  A.  12)  .  From  Figures  A.  10  -  A.  12 
and  7.3a  -  7.3c,  note  the  fundamental  frequency  for  the  full- 
scale  CRF  is  6.9  Hz;  whereas,  for  the  model,  the  fundamental 
frequency  is  6.6  Hz.  Again,  the  frequency  data  for  the  full- 
scale  facility  is  verified  through  experimental  results. 

At  this  point,  some  attention  should  be  given  to  the  low 
frequency  (f  «  1  Hz)  response  of  the  CRF  and  model.  Unlike 
the  blocked  case.  Figure  7.3  (and  A.  10  as  well)  indicates  that 
the  gains  are  not  zero  when  the  frequency  is  zero.  This  is 
simply  a  consequence  of  the  scale  of  the  plots.  Studies  were 
conducted  to  determine  if  the  analysis  does  go  to  zero  at  low 
frequency.  It  was  found  that  in  all  cases  considered,  the 
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gains  drop  very  rapidly,  on  the  order  of  10  dB/.OOl  Hz,  for 
frequencies  much  less  than  1. 

As  with  the  full-scale  results  in  Chapter  VI,  these 
results  show  that  the  peaks  reduce  with  increasing  end  radius. 
The  data  show  good  agreement  with  theory  although  the 
agreement  is  not  as  good  as  with  the  blocked  line  case. 

7.2.3  Screens  and  Honeycomb:  R-  -  1/2R3.  Figures  7.4a 
through  7.4c  illustrate  the  model  response  under  the  following 
conditions.  The  model  is  configured  for  basic  length,  screens 
and  honeycomb  for  the  flow  conditioning  elements,  and  an  end 
radius  of  ^3.  The  mass  flow  is  due  only  to  the  sinusoidal 
input  signal  and  is  computed  as  2.3  X  10~‘  slug/s.  The 
corresponding  full-scale  CRF  plots  are  located  in  Appendix  A 
in  Figures  A. 13  through  A. 15. 

Figures  7.4a  shows  the  gain  data  again  falling  relatively 
close  to  the  curves  below  30  Hz,  although  when  compared  with 
the  same  curve  for  the  blocked  and  MR3  cases,  the  data  and 
analyses  do  not  agree  as  closely.  Again,  this  is  a  direct 
result  of  the  end  impedance  equation  used  for  the  model  and 
indicates  the  end  impedance  equation  (2-15)  is  less  accurate 
for  the  more  open  lines. 

By  increasing  the  radius  to  ^3  from  ’4R3  in  the 
analytical  solution  for  the  model,  the  fundamental  peak  drops 
from  near  20  dB  (Figure  7.3a)  to  just  above  0  dB  (Figure 
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7.4a).  The  other  transducer  locations^  however/  are  not 
affected  as  greatly  by  Increasing  the  end  radius.  For  PR3/PS/ 
the  fundamental  peak  Increases  5  dB  from  the  ^3  case  (Figure 
7.3b)  to  the  Vsr3  case  (Figure  7.4b)  and  the  fundamental  peak 
for  PR2/PS  rises  5  dB  from  the  MR3  case  (Figure  7.3c)  to  the 
^3  case  (Figure  7.4c).  Primarily/  these  results  show  that 
the  end  impedance  dependence  of  the  line  is  a  function  of 
location.  In  particular/  the  end  impedance  strongly  affects 
the  end  transducer  but  does  not  have  a  big  effect  on  the 
upstream  transducers.  This  is  primarily  due  to  the  fact  that 
the  end  impedance  for  the  upstream  transducers  is  more 
strongly  dependent  on  the  impedance  mismatch  at  the  junction 
of  sections  2  and  3  of  the  model  than  on  the  orifice  radius. 
This  impedance  at  the  junction  of  sections  2  and  3  is  high  due 
to  the  smill  radius  of  section  3  as  compared  to  the  radius  of 
section  2.  In  effect,  the  end  impedance  for  sections  1  and  2 
is  nearly  that  of  a  blocked  line.  Thus,  changing  the  end 
radius  has  a  minimal  effect  on  the  signal  at  the  upstream 
transducer  locations. 

7.2.4  Screens  and  Honeycomb:  Open  Line  Response.  Figures 
7.5a  through  7.5c  illustrate  the  fully  open  behavior  of  the 
model.  The  full-scale  plots  are  located  in  Appendix  A  as 
Figures  A. 16  through  A.  18.  This  set  of  results  is  for  the 
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3/16  in,  m  =  2.1  X  10‘®  slug/s,  basic  length,  screens  and  honeycomb. 


same  configuration  as  the  last  three  sets  of  results,  except 
the  end  radius  Is  now  fully  open  and  equal  to  R3. 

Kote  first  that  Increasing  the  end  radius  to  R3  causes 
the  data  points  to  fall  closer  to  the  curve  than  in  the  two 
orifice-terminated  cases.  In  fact,  the  data  fall  on  the  curve 
about  as  well  as  for  the  blocked  case.  This  is  due  to  the 
fact  that  the  open  line  end  impedance  is  zero;  thus,  no 
approximations  are  required. 

Comparing  Figures  7.5a  and  7.4a,  the  fundamental  peak 
drops  10  dB  from  the  V4R3  case  to  the  fully  open  case. 
However,  the  gains  for  PR3  (Figure  7.4b  and  7.5b)  and  PR2 
(Figures  7.4c  and  7.5c)  change  very  little,  indicating  these 
two  locations  are  relatively  unaffected  by  radius  increases 
beyond  ^3 . 

The  full-scale  results  are  similar  to  the  model  results 
for  PR2  and  PR3  in  that  they  change  very  little  when  the 
radius  is  increased  from  ^3  to  R3.  Also,  for  the  full-scale 
analysis,  PR4/PS  changes  very  little  as  well.  This  appears  to 
contradict  the  results  of  Chapter  VI  where  it  was  noted 
several  times  that  the  gain  of  the  line  overall  facility  is 
highly  dependent  on  the  end  radius.  The  reasons  for  this 
seeming  discrepancy  are  explained  i^i  the  next  section. 

7.2.5  End  Conditions.  Transducer  Locations,  and  Gain. 
First,  it  is  noted  that  PR4/PS  is  not  the  gain  of  the  whole 
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system  (PEND/PIN)  .  PR4  is  not  located  exactly  at  the  exit  and 
PS  is  not  precisely  located  at  the  inlet.  It  turns  out  that, 
for  the  full-scale  facility,  these  locations  are  very 
important  to  the  gains  when  the  end  radius  is  greater  than 
zero. 

Second,  the  mass  flows  in  the  model  are  so  low  that  the 
resistance  at  the  exit  is  nearly  zero  for  all  cases  terminated 
with  an  orifice.  Recall  that  the  end  impedance  equation  (2- 
15)  is  directly  proportional  to  mass  flow  and  inversely 
proportional  to  area.  For  the  full-scale  analysis,  the  mass 
flow  is  v^ry  low  and  the  area  is  very  large.  Under  these 
conditions,  PR4/PS  is  far  enough  upstream  of  the  exit  to  be 
unaffected  by  the  radius  change  at  the  exit  and  therefore  the 
response,  under  all  of  the  orifice-terminated  conditions,  is 
similar  to  the  open  line  response. 

In  contrast,  PEND/PIN  is  significantly  affected  by  the 
radius  changes  as  this  location  is  right  at  the  exit.  Figures 
A. 19  through  A. 21  show  the  full-scale  CRF  PEND/PIN  response 
for  the  MR3,  ^3  and  fully  open  cases  presented  in  this 
chapter.  The  end  impedances  are  extremely  low  for  these  cases 
(3  X  10-’  Ibf-s/ft^  for  the  ^R3  case,  2  X  lO'®  Ibf -s/ft®  for  the 
^3  case,  and  1  X  10"’  Ibf  -s/ft®  for  the  R3  case)  . 

This  information  is  important  and  will  be  further 
discussed  in  section  7.5.1.  For  now;  however,  the  discussion 
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will  return  to  the  lab  model.  In  the  next  section,  the  second 
group  of  no  throughflow  results  are  presented. 

7.3  Model  Response  with  Washers.  Three  sets  of  figures, 
each  representing  a  different  end  radius,  are  presented  in 
this  section  for  the  basic  length  configuration  with  washers 
installed  for  the  flow  conditioning  elements.  Again, 
throughflow  is  provided  only  by  the  input  signal.  For  all 
three  figure  sets,  the  mass  flow  is  2.3  X  10~‘  slug/s. 

The  purpose  of  this  analysis  is  to  determine  if  the  area 
reduction  technique,  used  in  the  analysis  to  model  the 
impedance  of  the  flow  conditioning  elements,  correctly 
predicts  the  gains  and  phase  shift  at  the  three  transducer 
locations  in  the  laboratory  model. 

Figure  7 . 6  shows  the  model  response  at  the  end  transducer 
(PR4)  for  washers  installed  in  section  2  with  the  end  radius 
open  to  MR3.  This  plot  can  be  compared  with  the  identical 
screen  and  honeycomb  results  in  Figure  7.3a.  The  data  for  the 
screens  and  honeycomb  configuration  fall  closer  to  the 
analytical  curve  for  PR4/PS  than  the  data  for  the  washer 
configuration  at  the  same  location.  The  other  two  transducer 
locations  show  the  same  good  agreement  with  theory  for  both 
the  washers  and  the  screens  and  honeycomb  (J^pendix  A,  Figures 
A. 22  through  A. 23)  . 
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Figure  7.6.  Model  Frequency  Response  (PR4)  with  Washers  for  R^  =  ViR3 .  Conditions: 
=  3/64  in,  ih  =  2.3  X  10‘*  slug/s,  basic  length,  washers. 
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Figure  7.7  shows  the  model  response  for  washers  in 
section  2  and  the  end  radius  open  to  ^iR3.  This  plot  is 
compared  to  the  identical  screens  and  honeycomb  results  in 
Figure  7.4a.  For  this  configuration,  the  washer  results  show 
better  agreement  with  theory  for  PR4/PS  than  the  screens  and 
honeycomb  results.  Near  the  flow  conditioners  (PR3  and  PR2) 
both  configurations  again  show  excellent  agreement  with  theory 
(Figures  A. 24  and  A. 25) . 

Finally,  Figure  7.8  illustrates  the  model  response  for 
PR4  with  washers  in  section  2  and  the  end  radius  fully  open  to 
R3.  These  plots  compare  with  Figure  7.5a.  In  this  case, 
neither  configuration  shows  better  agreement  with  theory  than 
the  other  for  any  of  the  three  transducer  locations.  PR2/PS 
and  PR3/PS  plots  are  in  Appendix  A,  Figures  A. 26  and  A. 27. 

The  data  appear  to  be  independent  of  the  flow 
conditioning  elements  as  no  one  configuration  is  better 
matched  to  theory  than  another.  The  reason  the  flow 
conditioning  elements  do  not  affect  the  solution  is  because 
the  impedance  presented  by  these  elements  is  small  in 
comparison  to  the  area  changes  that  occur  between  sections  of 
the  facility.  Under  no  throughflow  conditions,  since  the  area 
reductions  are  small,  the  flow  straightening  elements  are 
virtually  nonexistent  and  do  not  affect  the  analytical  model. 
Unfortunately,  this  conclusion  sheds  no  light  on  whether  or 
not  the  area  reductions  are  correct  models  for  the  flow 


7-28 


Model  Frequency  Response  (PR4)  with  Washers  for  R^  =  R3 .  Conditions 
=  2.3  X  10'®  slug/s,  basic  length,  washers. 


conditioning  elements.  Data  was  taken  with  the  screen  only 
flow  conditioning  element  section;  however,  the  results  differ 
little  from  the  screen  and  honeycomb  results.  Therefore,  the 
screen  plots  are  not  presented. 

_ Model  Response  with  Throuahflow 

In  this  section,  model  results  are  presented  for  the 
basic  length  configuration  with  screens  and  honeyccmib 
installed  for  the  flow  conditioning  elements  and  throughflow 
present  in  the  system.  Three  experimental  runs  are  presented, 
each  for  a  different  end  radius.  The  experimental  conditions 
are  similar  to  the  conditions  for  the  results  shown  in  section 
7.2.2  -  7.2.4  except  for  the  additional  variedsle  of 
throughflow.  For  the  results  presented  in  the  following 
paragraphs,  signal  propagation  is  affected  by  the  velocity  of 
the  throughflow.  The  velocity  of  the  waves  propagating  in  the 
direction  of  the  throughflow  are  traveling  at  the  speed  of 
sound  plus  the  throughflow  velocity.  Waves  propagating 
against  the  direction  of  throughflow  are  traveling  at  a 
velocity  of  the  speed  of  sound  minus  the  throughflow  velocity. 

The  results  in  Figures  7.9a  through  7.9c  are  for  an  end 
radius  of  ^3  and  a  mass  flow  rate  of  8.8  X  10~’  slug/s.  This 
flow  rate  corresponds  to  a  velocity  of  31.7  ft/s  in  section  3 
of  the  model.  This  velocity  corresponds  to  a  mass  flow  of  28 
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l]:xn/s  in  the  full-scale  facility  and  is  therefore  similar  to 
the  20  Itmi/s  analyses  presented  in  Chapter  VI. 

Note  the  fairly  good  agreement  of  theory  with  the  data 
for  PR4  and  PR2.  PR3  data;  however,  does  not  agree  with 
theory  well  at  all.  Recalling  that  PR3  is  located  just  aft  of 
the  flow  conditioning  elements,  it  is  probed>le  that  the  flow 
is  disturbed  at  this  location  and  the  data  not  very  reliable. 
Since  the  transducers  are  configured  to  read  any  pressure 
changes  in  the  line,  any  fluctuating  pressure  signal  will  be 
picked  up  by  PR3.  It  is  well  known  that  screens  cause  a  large 
pressure  drop  (Rae  and  Pope:74)  .  It  is  likely  that  the 
pressure  drop  across  the  flow  conditioning  elements  is 
unsteady  and  this  fluctuation  is  picked  up  by  the  transducer. 
With  fluctuations  from  the  screen  and  throughflow 
interactions,  the  input  signal  is  harder  to  distinguish  when 
input  into  the  spectral  analysis  program.  This  is  primarily 
a  local  effect  since  PR2  data  agrees  well  with  theory  and  PR4 
data  is  more  dependent  on  end  condition  than  flow  quality. 

Comparing  Figure  7 . 9a  with  the  results  from  the  same 
configuration  without  throughflow.  Figure  7.3a,  it  is  noted 
that  for  PR4/r5?,  the  data  for  throughflow  agrees  with  theory 
about  as  well  as  the  data  for  the  case  without  throughflow. 
This  is  also  true  for  PR2  as  is  clear  when  comparing  Figures 
7.8c  and  7.3c.  Also,  the  fundamental  frequency  is  reduced 
from  6.6  Hz,  for  the  no  throughflow  case,  to  6.3  Hz  for  the 
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results  with  throughflow.  Further,  the  peak  levels  are  higher 
with  throughflow  than  without.  Both  the  lower  frequency  and 
the  higher  peaks  are  indicative  of  higher  end  impedance.  This 
is  indeed  the  case,  as  the  end  impedance  for  the  throughflow 
case,  computed  from  Equations  (2-15)  and  (2-16) ,  is  500  times 
higher  than  for  the  no  throughflow  case.  The  mass  flow,  then, 
does  affect  both  the  peaks  and  frequencies  of  the  physical 
model . 

Figures  7.10a  through  7.10c  show  the  data  and  analysis 
for  the  ^3  case.  Again  the  data  shows  good  agreement  with 
theory  and  when  compared  to  the  identical  no  throughflow  cases 
in  Figures  7.4a  through  7.4c.  The  analysis  for  PR4/PS  are 
slightly  different  for  the  no  throughflow  case  in  Figure  7.4a 
than  for  the  throughflow  case  of  7.10a.  Specifically,  the 
throughflow  PR4/PS  curves  illustrate  higher  gains.  The  fact 
that  the  data  agrees  with  theory  in  both  throughflow  and  no 
throughflow  cases  means  that  the  model  is  sensitive  to  the 
changes  in  conditions  of  the  flow.  Note  that  in  the 
throughflow  case,  PR4/PS  data  matches  more  closely  with  the 
analysis  than  in  the  case  without  throughflow.  Also,  PR3  and 
PR2  data  more  closely  agree  with  theory  than  in  the  ^R3  case, 
even  though  the  theoretical  curves  from  the  two  locations 
change  very  little  from  the  no  throughflow  curves. 

For  this  experiment,  the  mass  flow  is  1.73  X  10“^  slug/s 
and  the  velocities  in  sections  1,  2,  and  3  are  2.4  ft/s,  9.7 


7-36 


3/32  in,  ih  =  1.73  X  10'^  slug/s,  basic  length,  screens  and  honeycomb. 
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Figure  7.10c.  Model  Throughflow  Frequency  Response  (PR2)  for  R^  =  V4R3 .  Conditions: 
Rc  =  3/32  in,  ih  =  1.73  X  10‘*  slug/s,  basic  length,  screens  and  honeycomb. 


ft/a,  and  69.2  ft/s,  respectively.  The  same  section  3 
velocity  in  the  CRF  corresponds  to  a  mass  flow  of  60  Ibm/s  in 
the  full-scale  facility. 

Compared  with  the  conditions  for  the  results  in  Figures 
9a  through  9c,  the  mass  flow  is  increased  and  the  end 
condition  more  open.  This  lowers  the  real  part  of  the  end 
impedance  [Equation  (2-15) ] ;  however,  it  is  still  680  times 
higher  than  the  end  impedance  for  the  case  without  flow.  As 
was  the  case  with  the  ^3  end  radius  results,  the  fund2uiental 
frequency  drops  when  flow  is  introduced.  In  this  case,  the 
fundamental  frequency  is  6.4  Hz  for  the  results  with  flow  and 
6.7  Hz  for  the  results  without  flow. 

The  last  three  Figures,  7.11a  through  7.11c,  are  for  an 
end  radius  of  ^3.  In  this  case,  the  mass  flow  is  2.2  X  10'* 
slug/s  and  the  velocity  in  section  3  is  90.5  ft/s.  This 
velocity  is  most  representative  of  flow  conditions  that  would 
be  seen  in  the  full-scale  facility.  Specifically,  this 
velocity  in  section  3  corresponds  to  a  80  Ibm/s  flow  rate  in 
the  full-scale  facility. 

These  results  show  that  as  the  end  orifice  is  increased 
from  ^3,  the  analytical  and  experimental  gains  decrease  for 
PR4/PS.  The  fundamental  frequency  increases  from  6.4  Hz,  for 
the  ^3  case,  to  6.7  Hz  for  the  ^3  case. 

In  general,  the  fundamental  frequency  data  in  Figures  7.9 
through  7.11  do  not  agree  as  well  with  theory  as  the  no 
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throughflow  cases.  Specifically,  the  end  transducer  data  do 
not  fall  on  the  curve  well  at  the  fundamental  frequencies.  At 
low  throughflow  rates  and  high  impedance  (Figure  7.9),  the 
theory  overestimates  the  fundamental  peak.  At  higher 
throughflow  rates,  the  theory  underestimates  the  fundamental 
peak.  This  is  a  result  of  the  analysis  not  fully  representing 
the  physical  system.  Losses  due  to  friction  are  computed  but 
may  not  be  exact,  resulting  in  a  mismatch  between  theory  and 
the  physical  model  data  that  is  amplified  when  throughflow  is 
introduced.  Also,  the  end  impedance  model  is  only  a  best 
attempt  at  representing  the  line  end  and  is  most  likely  not 
the  true  representation  of  the  end  condition. 

The  frequency  response  results  presented  for  the  model 
support  the  analysis  of  the  full-scale  facility  and 
illustrates  some  of  the  points  made  in  Chapter  VI  concerning 
gain  and  frequency  dependence  on  end  impedance.  Most 
importantly,  the  laboratory  experiments  verified  that  the 
fundamental  frequency  is  near  6.5  Hz  for  the  basic  length  case 
with  a  wide  variety  of  flow  conditions  and  end  configurations. 

Also,  the  model  verified  some  trends  in  the  analyses 
presented  in  Chapter  VI.  Specifically,  under  no  throughflow 
conditions,  opening  up  the  end  radius  decreased  the  gains  and 
increased  the  fundamental  frequency  for  the  model  in  the  same 
manner  as  the  full-scale  facility.  This  was  demonstrated  both 
in  the  analysis  and  in  the  data.  As  with  the  full-scale 
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facility,  the  gains  increased  with  increasing  throughflow  and 
the  fundamental  frequencies  increased  with  increasing  end 
radius . 

7.5  Transient  Response  of  the  Model 

The  transient  response  tests  were  accoBg)lished  to 
determine  the  damping  characteristics  of  the  model.  For  the 
transient  response  tests,  the  model  configuration  includes 
screens  and  honeycomb,  basic  length,  and  em  end  radius  of  UR3. 
A  solenoid  valve  was  used  to  provide  a  positive  pressure  step. 

The  results  are  presented  in  four  figures.  Figure  7.12 
illustrates  the  long  time  response  for  a  0.08  psig  inpi^  step. 
For  this  pressure  input  and  end  condition,  the  model  reaches 
a  steady  state  condition  in  4.5  seconds.  For  larger  inputs, 
the  model  takes  longer  to  reach  a  steady  state  pressure. 

Figure  7.13  shows  the  unsteady  response  to  the  input 
signal  at  the  end  transducer,  PR4.  The  valve  opens  at  time 
t  ■■  -  0.06  seconds.  At  t  -  0  seconds,  the  pressure  is  felt  at 
the  end  transducer. 

Initially,  there  is  no  flow  in  the  model.  At  t  -  0,  the 
step  function  from  the  solenoid  valve  reaches  the  end  of  the 
line  and  PR4  senses  the  pressure  rise.  Immediately  after  the 
pressure  front  reaches  the  end  of  the  line, ' the  sudden  flow 
constriction  of  the  orifice  causes  a  reflection  of  the  front 
and  the  pressure  immediately  drops.  The  reflected  wave 
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Figure  7.12:  Long-Time  Model  Response  (PR4)  to  a  Positive  Step 
Input.  Conditions:  R^  *  ^<R3,  basic  length,  screens  and 
honeycomb . 
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Figure  7.13:  Model  Transient  Response  (PR4)  to  a 
Small, Positive  Step  Input.  Conditions:  Rc  *  ^4^,  basic 
length,  screens  and  honeycomb. 


returns  to  the  inlet,  slightly  damped  £r<»B  viscous 
interactions  with  the  boundary  layer,  and  reflects  back  to  the 
end  where  PR4  again  picks  up  its  presence  at  time  0.12 
seconds.  These  reflections  continue  until  the  wave  is  fully 
damped  at  t  ■■  0.7  seconds.  The  pressure  fluctuations  from  the 
wave  reflecting  back  and  forth  between  the  inlet  and  exit  are 
termed  "ringing" . 

This  pressure  step,  of  0.08  psig,  is  on  the  order  of  the 
pressure  level  of  the  signal  introduced  by  the  pneumatic 
signal  generator  for  the  frequency  plots  presented  in  the 
previous  two  sections.  Figure  7.13  shows  that  for  the  small 
signal  inputs  in  the  frequency  d<xaain,  the  associated  ringing 
is  not  very  significant.  In  fact,  the  steady  state  pressure 
is  not  yet  reached  by  the  time  the  ringing  is  damped  out. 
This  means  that,  for  the  model,  the  transient  response  is 
controlled  by  the  time  required  to  reach  steady  state 
pressures  and  not  by  the  ringing.  This  also  means  the  model 
is  well  damped. 

Figures  7.14  and  7.15  illustrate  that  the  model  also 
exhibits  damped  behavior  for  larger  pressure  steps.  Figure 
7.14  shows  the  ringing  associated  with  a  pressure  step  larger 
than  0.45  psig.  Again,  the  ringing  ceases  in  a  short  time; 
here,  within  0.6  seconds.  As  with  the  small  step,  the  ringing 
damps  out  before  the  steady  state  pressure  is  reached. 
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Figure  7.14:  Model  Transient  Response  (PR4)  to  a  Large 
Positive  Step  Input.  Conditions:  Re  -  ^<R3,  basic  length, 
screens  and  honeycomb,  AP  >  .45  psig. 
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Figure  7.15:  Model  Transient  Response  (PR4)  to  a  Very  Large 
Positive  Step  Input.  Conditions:  Re  ~  basic  length, 
screens  and  honeycomb,  AP  >  1  psig. 


Figure  7.15  illustrates  the  ringing  for  an  even  larger 
pressure  step,  greater  than  1  psig.  Again  the  ringing  ceases 
within  0.6  seconds  and  before  the  steady  state  pressure  is 
achieved . 

These  results  show  that  the  model  is  well  damped  and  the 
signal  size  may  affect  the  ringing  time.  In  Figure  7.13,  with 
a  small  step,  the  ringing  damps  out  in  0.7  seconds;  whereas, 
with  the  larger  pressure  steps  in  Figures  7.14  and  7.15,  the 
ringing  damps  out  in  0.6  seconds.  There  is  not  enough 
evidence  to  make  any  broad  statements  about  damping  levels  and 
signal  sizes,  but  it  is  noted  that  the  results  tend  to  show 
the  signal  size  is  of  some  importance;  particularity,  the 
larger  signals  damp  out  more  rapidly  than  the  smaller  ones. 

Negative  pressure  step  experiments  were  also  run; 
however,  the  data  is  not  presented  here.  Essentially,  no 
ringing  was  observed  for  the  negative  step  and  therefore,  the 
data  provide  no  additional  insight  into  the  damping  of  signals 
in  the  model. 

The  transient  analysis  on  the  model  can  be  loosely 
associated  with  the  full-scale  facility  damping 
characteristics.  In  general,  the  model  is  well  damped  and, 
because  of  this,  one  might  expect  the  full-scale  facility  also 
be  well  damped.  Bear  in  mind;  however,  the  radii  of  the  model 
sections  are  120  times  smaller  than  the  radii  of  the  full- 
scale  facility.  Boundary  layer  interactions  in  the  model  will 
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be  much  more  significant  than  in  the  full-scale  facility/  due 
to  the  smaller  radius.  As  such,  reflections  will  be  damped 
out  faster  in  the  model  than  in  the  full-scale  facility.  In 
any  case,  it  is  not  unreasonable  to  predict,  on  the  basis  of 
the  model  transient  results,  that  the  full-scale  CRF  damping 
is  high.  Further,  it  is  reasonable  to  expect  the  input  signal 
size  is  a  factor  in  damping  characteristics  of  the  full-scale 
facility,  ju£t  as  it  seems  to  be  in  the  characteristics  of  the 
model  damping. 

As  previously  mentioned,  the  CRF  damping  is  of 
significant  importance  to  the  control  system  design.  At  this 
time,  the  discussion  postponed  in  Chapter  VI  is  picked  up 
again . 

7  ^6  CRF  Damping 

Attempts  were  made  by  the  controls  designer  to  correlate 
the  CRF  frequency  response  curves  presented  in  Chapter  VI  with 
the  system  damping  (Hull : 1993) .  Unfortunately,  the  attempt 
was  unsuccessful.  Hull's  time  response  computations,  based  on 
the  full-scale  frequency  response  results,  yielded  a  very 
lightly  damped  system.  Several  years  of  experience  running 
the  CRF  has  proven  the  system  to  be  well  damped  (Hull:  1993) 
and  consequently,  indicate  error  in  the  amplitude  results  of 
this  analysis.  Several  sources  of  error  are  the  end  impedance 
model,  the  velocity  adjustments  to  the  propagation  operator 
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and  characteristic  impedance,  and  the  small  pressure  signal 
assumption.  First,  the  errors  associated  with  the  end 
Impedance  model  are  discussed. 

7.6.1  End  Impedance  Model  Errors.  Reviewing  the  results 
of  the  experimental  runs,  it  is  clear  that  the  analysis  and 
experimental  data  match  up  quite  well  for  the  low  flow  rate 
runs.  However,  it  is  also  clear  that  the  analysis  does  not 
correctly  predict  gains  when  the  mass  flow  is  increased. 
Since  the  experimental  data  yielded  good  agreement  with  theory 
at  known  end  impedance  conditions  (blocked  and  fully  open) , 
this  indicates  that  the  analysis  is  lacking  some  element 
connected  to  the  mass  flow  conditions. 

The  gains  are  dependent  on  mass  flow  in  many  ways;  the 
most  obvious  of  which  is  the  end  impedance.  Recalling  the 
results  from  this  Chapter  and  Chapter  VI,  the  end  impedance 
has  a  significant  effect  on  the  amplitudes  of  the  peaks.  The 
overall  gain  of  the  CRF  is  highly  dependent  on  the  location  of 
the  end  of  the  line  as  well  as  on  the  impedance  model. 

As  was  seen  in  section  7.2,  the  location  of  the  end  of 
the  line  can  have  a  significant  affect  on  the  gains.  In 
section  7.2,  it  was  shown  that  under  low  flow  conditions, 
PR4/PS  gains  are  significantly  different  from  PEND/PIN  gains. 
PEN  is  located  just  1  7/8  in  downstream  of  PR4  and  PIN  is  just 
0.333  in  upstream  of  PS.  More  study  could  be  accomplished  to 
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better  define  the  end  and  inlet  locations  of  the  CRF,  possibly 
resulting  in  different  gain  curves. 

The  end  impedance  model  used  for  the  full-scale  analysis 
was  determined  from  theory  and  is  only  a  best  attempt  at  a 
model.  The  solution  was  simplified  by  ignoring  the  compressor 
in  all  ways  except  size.  Further  study  could  be  performed  to 
more  completely  represent  the  compressor  in  the  end  impedance 
model  and  predict  the  peak  amplitudes  more  accurately.  If 
this  study  is  performed/  the  velocity  effects  on  the 
propagation  operator  and  characteristic  impedance  should  also 
be  refined  as  discussed  in  the  next  section. 

7.6.2  Errors  with  Velocity  Effects.  The  equations  used 
to  include  velocity  effects  [(2-17)/  (2-18)/  and  (2-19)]  are 
only  an  approximation.  In  this  analysis,  the  approximation 
was  used  to  simplify  the  calculations  and  consequently,  may 
have  contributed  to  the  problem  of  high  gains.  Katz,  Hausner, 
and  Eisenberg  (1974:282-285)  present  a  more  precise  solution 
for  including  high  velocities  that  requires  simultaneous 
solution  of  the  transmission  line  equations. 

Including  the  velocity  effects  using  the  more  precise 
solution  may  help  reduce  the  gains.  However,  there  is  no  way 
to  estimate  the  extent  to  which  the  gains  will  decrease,  if  at 
all. 
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order  to  truly  present  a  transfer  function  useful  to  the 
controls  designer,  the  governing  equations  must  be  refined  to 
include  large  amplitude  signals.  The  controls  designer  needs 
to  input  large  pressure  changes,  of  magnitude  2-5  psig.  The 
analysis  presented  in  this  paper  assumes  small  signals  of 
order  2  or  3  lower  than  the  required  signal  size. 

The  results  of  the  transient  tests  with  the  model 
indicate  that  large  pressure  signals  may  damp  out  faster  than 
small  signals.  Kantola  (1971:27^-281)  has  also  shown  this  to 
be  true.  As  such,  an  analysis  for  large  input  signals  should 
result  in  lower  gains  than  with  this  small  signal  analysis. 
In  order  to  include  this  effect  in  the  analysis;  however,  the 
Navier-Stokes  equations  cannot  be  linearized. 

Another  assumption  that  is  most  likely  violated  is  the 
plane  wave  assumption.  In  lines  with  radii  small  compared  to 
the  length,  such  as  in  the  model,  this  assumption  is  valid. 
However,  in  the  large  radii  of  the  CRF,  the  signals  traveling 
the  69.8  ft  length  of  the  line  can  develop  disturbances  in  the 
radial  direction,  resulting  in  waves  that  are  more 
spherically-shaped  than  plane-shaped. 


The  primary  objective  of  this  effort  was  to 
determine  the  fundamental  frequency  of  the  CRF  and  this  was 
successfully  acconplished.  The  resonant  frequencies 
calculated  for  the  full-scale  facility  were  checked  against 
basic  pipe  resonance  equations  and  verified  through 
experimental  results. 

8.2  Conclusion 

All  of  the  analyses  place  the  fundamental  frequency 
between  5.5  Hz  and  6.9  Hz.  The  peak  amplitudes  for  the  full 
scale  facility  are  questionable;  however,  the  focus  of  this 
effort  was  to  find  the  resonant  frequencies  and  that  has  been 
acconplished. 

As  such,  it  is  recommended  the  control  system  be 
designed,  if  practical,  to  operate  at  a  frequency  well  below 
5.5  Hz.  This  frequency  occurs  at  the  worst  case  conditions 
(Run  10.1,  extended  length,  blocked  line). 

8.3  Recommendations 

Future  studi.es  on  the  CRF,  to  determine  the  true  gains, 
are  possible  if  the  effort  is  required  for  the  control  system 
design.  Recommendations  on  iirproving  the  analysis  for  the  CRF 
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were  presented  in  Chapter  VII  and  could  be  used  to  continue 
this  study.  At  present;  however,  the  control  system  design  is 
proceeding  without  the  correct  gain  information  and,  at  this 
point,  it  appears  that  future  studies  on  the  CRF  dynamic 
response  are  not  required. 

However,  some  basic  studies  could  be  accomplished  to 
broaden  the  knowledge  base  of  fluid  transmission  line  theory. 
A  simple  fluid  transmission  line  model  could  be  constructed, 
from  a  single  line  of  constant  diameter,  to  study  the  effects 
of  end  impedance  or  the  iiipedance  presented  by  flow 
conditioning  elements.  In  this  manner,  the  geometry  variables 
can  be  reduced,  from  those  in  the  CRF  model,  so  the  research 
can  concentrate  on  modeling  one  physical  property  of  the  fluid 
transmission  line. 

The  simple  model  could  be  constructed  from  rigid  tubing 
and  have  transducers  located  at  the  inlet  and  exit  of  the 
line.  The  frequency  response  of  the  simple  model  could  be 
determined  from  a  variety  of  input  signals  and  flow 
conditions.  From  these  tests,  more  baseline  data  would  be 
gathered  to  use  in  future  frequency  response  analyses  of  the 
type  presented  in  this  thesis. 
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The  following  pages  contain  additional  data  supporting 
Chapters  VI  and  VII.  Pages  A-2  through  A-10  contain  data  in 
support  of  Chapter  VI.  Pages  A-11  through  A-28  support 
Chapter  VII . 


Run  8.6 


Response  for  Basic  Length  with  20  Ibm/s  Mass  Flow  Rate.  Conditions 
825  ft  =  .  Fundamental  Peak  at  6.9  Hz. 


Run  10.6 


Run  9.7 


Run  10.7 


Run  8.9 


Figure  A. 7:  CRF  Response  for  Basic  Length  with  250  Ibm/s  Mass  Flow  Rate.  Conditions: 
Run  8.9,  K  =  0.825  ft  =  ViR3.  Fundamental  Peak  at  6.9  Hz. 
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Run  10.9 


8 


S 


8 


8 


o 


N 

I 


f 

u. 


4) 

9c- 


in 


b 

“  ij 

n  S 

cu 


I 


^  <0 

o  d 

in  S 

nj  c 


£l 

»£ 

.C 

a  • 

cS 
-o  " 
0) 


I 


in 

r» 

00 


U'-' 
®  II 

(uaT 

(Q 

d  «• 
O  oi 
O.  • 
IQ  O 
0)  <-l 
OC 

d 

bn  d 
Qi:  a: 

o 


a\  to 

<  o 

'H 
<U  4J 

Vi  rl 

3  TJ 
O)  d 
-H  O 
bt  U 


RunEISFS 


•H  • 
O 


gure  A.  11:  CRF  No  Throughflow  Response  (PR3)  for  R^  =  HR3 .  Conditions: 
4375  ft,  ih  =  2.2  X  10'®  slug/s,  basic  length. 
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Throughflow  Response  (PR2)  for  Rc  =  HR3 .  Conditions 
;  10‘*  slug/s,  basic  length. 
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Figure  A.  15:  CRF  No  Throughflow  Frequency  Response  (PR2)  for  R^  =  WR3  .  Conditions: 
R,  =  0.825  ft,  A  =  2.3  X  10'*  slug/s,  basic  length. 


Figure  A.  16:  CRF  No  Throughflo' 
R.  =  1.65  ft,  ih  =  2.1  X  10-*  sl\ 


iroughflow  Frequency  Response  (PR3)  for  R^  =  R3 .  Conditions: 
X  10'*  slug/s,  basic  length. 
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Throughflow  Frequency  Response  (PR2)  for  R^  *  R3 .  Conditions 
10'®  slug/s,  basic  length. 
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ure  A.  19:  CRF  No  Throughflow  Response  (PEND)  for  R,.  =  KR3 .  Conditions: 
375  ft,  ih  =  2.2  X  10'*  slug/s,  basic  length. 
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Throughflow  Response  (PEND)  for  Rc  =  Vfil3 .  Conditions 
10'®  slug/s,  basic  length. 


Figure  A. 21:  CRF  No  Throughflow  Frequency  Response  (PEND)  for  R^  =  R3 .  Conditions: 
R_  =  1.65ft,  ih  =  2.1X  10‘‘  slug/s,  basic  length. 
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Figure  A. 23;  Model  Frequency  Response  (PR2)  with  Washers  for  R^  =  MR3 .  Conditions: 
Rc  =  3/64  in,  ih  =  2.3  X  10'*  slug/s,  basic  length,  washers. 


length 


)  with  Washers  for  =  R3 .  Conditions 
length,  washers. 


This  appendix  presents  a  chronological  overview  of  the 


development  of  the  computer  program  used  to  solve  the 
Nichols  equations.  The  program  was  developed  first  for  a 
constant  diameter  line  with  a  blocked  end.  All  of  the 
geometry  and  flow  effects  were  added  in  increments  to 
determine  how  the  solution  is  affected  by  each. 

First,  the  Nichols  solutions  for  a  constant  diameter 
blocked  line  were  solved  and  verified.  Next,  a  study  on  the 
effect  of  scaling  the  radius  for  the  model  dimensions  was 
accomplished.  The  solution  was  then  modified  to  account  for 
cascaded  lines  of  different  diameter  and  orifice 
terminations.  Finally,  the  flow  conditioning  elements,  the 
pressure  equalization  pipe,  and  throughflow  were  added  to 
the  solution.  Some  of  the  above  developments  are  discussed 
here.  Others,  such  as  the  effect  of  end  radius  and 
throughflow,  are  covered  in  detail  in  Chapters  VI  and  VII. 

B.l  Constant  Diameter.  Blocked  Line  and  Radius  Scaling 

The  solution  for  a  constant  diameter  line  of  20  ft 
diameter  and  69.8  ft  in  length  is  shown  in  Figure  B.l.  The 
gains  smoothly  attenuate  through  the  frequency  range;  but, 
they  are  high  (70  dB  for  the  fundamental  peak) .  In  Figure 
B.2,  a  constant  diameter  line  of  2  in  and  69.8  ft  in  length 


0  20  40  eo  W  100  120 
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Figure  B.2:  Computer  Output  for  a  2  in  Diameter  Line. 
Conditions:  Constant  Diameter,  Blocked,  69.8  ft  Long. 
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is  shovm.  Notice  the  gains  are  much  lower  but  the  shape  of 
the  curve  is  similar  to  Figure  B.l.  The  diameter  of  the 
lines  in  Figures  B.l  and  B.2  represent  section  1  of  the  CRF 
and  model,  respectively.  The  length,  69.8  ft,  of  both  lines 
represents  the  overall  length  of  the  CRF  and  model. 

Several  other  diameter  lines  were  studied  with  the 
constant  diameter  solution.  It  was  found  that  the  gains  are 
directly  dependent  on  the  radius  of  the  line  as: 

Gain{D)  =  Gain{20  ft  dia)  ^  .  (B-1) 

For  the  20  ft  diameter  line,  the  fundamental  peak  gain  is 
2400  (70  dB)  and  the  2  in  line  gain  is  20  (26  dB) .  The 
gains  are  related  by  the  above  equations.  The  gains  of 
other  radius  lines  studied  also  follow  equation  B-1  as  can 
be  shown  from  the  information  in  Table  B.l. 

The  radius  scaling  also  affects  the  fundamental 
frequency  as  is  clear  from  Table  B.l.  The  model  scale  is 
1:120  for  sections  1  and  2  and  the  pressure  equalization 
pipe.  According  to  the  table,  the  model  fundamental 
frequencies  should  be  slightly  lower  than  the  full-scale 
CRF. 

The  fundamental  frequency  computed  for  the  full-scale 
dimensions  of  20  ft  diameter,  4.08  Hz,  is  the  equal  to  the 
fundamental  frequency  calculated  from  the  basic  pipe 
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resonance  equation  for  a  blocked  line  (2-1)  of  69.8  ft  in 
length.  This  indicates  that  the  20  ft  diameter  line  is 
essentially  lossless,  since  Equation  (2-1.)  does  not  account 
for  viscous  effects.  The  smaller  radius  lines  have  lower 
fundamental  frequencies  than  the  20  ft  diameter  line, 
indicating  the  viscous  interactions  present  in  the  smaller 
lines  have  a  greater  affect  the  signal  propagation  than  in 
the  large  diameter  line. 

B.2  Cascaded  Lines  of  Different  Diameter. 

The  addition  of  the  10  ft  and  3.3  ft  diameter  sections 
to  the  20  ft  constant  diameter,  blocked  line  creates  a 
cascade  and  alters  the  solution  for  the  line.  Figure  B.3 


Table  B.l. 
Line .  The 

Effects 
speed  of 

of  Radius  Scaling  on 
sound  in  all  cases  is 

a  69.8  ft  Blocked 
1130.5  ft/s. 

Radius 

R 

(ft) 

Scale 

Fundamental 

Peak  Amplitude 
( nondimens ional ) 

Fundamental 

Frequency 

(Hz) 

10 

1:1 

2400 

4.08 

0.5 

1:20 

120 

4.03 

0.25 

1:40 

60 

4.01 

0.125 

1:80 

30 

3.97 

0.10 

1:100 

24 

3.94 

0.083 

1:120 

20 

3.92 

B-4 


Fretpjency.  f.  Hz 


Figure  B.3:  CRF  as  a  Three-Line  Cascade.  Conditions: 
Blocked,  Basic  Length. 


shows  the  blocked  line  with  a  sirt^jle  three-line  cascade. 
This  is  the  most  fundamental  solution  for  the  CRF  where 
section  1  is  of  20  ft  diameter,  section  2  is  of  10  ft 
diameter,  and  section  3  is  of  3.3  ft  diameter.  Shown  in 
Figure  B.3  is  the  solution  for  the  basic  length  line. 

Introduction  of  the  different  diameter  lines  produces 
an  impedance  mismatch  at  the  junctions  of  the  individual 
lines.  This  causes  the  fundamental  frequency  to  increase, 
in  this  case,  to  6.5  Hz  from  4.08. 


Changing  the  length  configuration  from  basic  to 
intermediate  decreases  the  fundamental  frequency  by  0.5  Hz. 
Another  0.5  Hz  drop  is  evident  by  changing  the  length 
configuration  from  intermediate  to  extended. 

Introduction  of  the  cascaded  lines  does  not  change  the 
peak  anplitude  of  the  fundamental  frequency,  but  it  does 
affect  the  harmonic  peaks  because  the  cascade  causes  signals 
to  reflect  at  the  junctions  of  the  individual  lines.  These 
reflected  signals  interact  with  each  other  and  the  input 
signal,  causing  the  attenuation  to  change  from  that  of  the 
constant  diameter  line  (Figure  B.l). 

B.3  Flow  Conditioning  Elements 

The  addition  of  the  flow  conditioning  elements  is 
analytically  modeled  as  area  reductions  in  the  locations 
where  they  are  mounted.  This  increases  the  number  of 
different  diameter  lines,  changing  the  cascade  from  3  to  13 
lines  of  different  diameter. 

The  solution;  however,  is  unaffected  by  the  addition  of 
the  flow  conditioners.  However,  because  the  throughflow  had 
not  yet  been  modeled,  the  flow  conditioning  elements  were 
left  in  the  program  for  all  subsequent  runs. 
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The  addition  of  the  pressure  equalization  pipe  to  the 
solution  for  the  cascaded  lines  changes  the  gain  curve  only 


slightly.  Figures  B.4  and  B.5  illustrate  the  blocked 
cascaded  line  (with  flow  conditioning  elements)  for  the  CRF 
with  and  without  the  pressure  equalization  pipe. 

The  solution  with  the  branch  modeled  shows  small  peaks 
at  30  and  111  Hz.  These  peaks  are  the  only  difference 
between  the  solution  with  and  without  the  branch.  Again, 
due  to  the  fact  that  throughflow  has  not  yet  been  accounted 
for,  the  pressure  equalization  pipe  was  left  in  the  model 
for  all  subsequent  runs . 

B.5  Sximmarv 

The  most  significant  effect  on  the  solution  for  the  CRF 
is  the  cascaded  lines  of  different  diameter  created  from 
sections  1,  2,  and  3.  The  solution  is  highly  dependent  on 
the  length  configuration;  however,  it  is  not  affected  by  the 
addition  of  flow  straightening  elements.  The  solution  is 
slightly  changed  with  the  addition  of  the  pressure 
equalization  pipe;  but  this  effect  is  at  higher  frequencies 
and  does  not  change  the  fundamental  frequency  at  all. 


Gain,  Rand/Pin,  db 


Figure  B.4:  CRF  with  Flow  Conditioning  Elements  (no  Pressure 
Equalization  Pipe) .  Conditions:  Blocked,  Basic  Length. 
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Pneumatic  Signal  Generator:  Manufacturer:  Bendix.  Includes 

Pressure  Regulator,  S/N  ENY-A296,  0  -  100  psig 

Lambda  Power  Supplies,  S/N  ENY-A295/6,  0  -  600  V,  0  -  200  mA 

Bendix  Voltage  Regulator,  S/N  ENY-A293 

Hewlett-Packard  Artplifier,  Model  450A,  S/N  ENY-A297 

Data  Acquisition  System: 

Nicolet  Multipro  Pedestal,  S/N  IBK9300114.  Includes  four 

Nicolet  Multipro  120  digitizer  cards,  S/N  IBM9300191/2/3/4 

House  Air  Mean  Pressure  Transducer 

Transducer:  Endevco  8510-B100,  S/N  20LG,  100  psig 

Signal  Conditioner:  Endevco  4423,  S/N  BA88 

Power  Supply:  Endevco  4425,  S/N  AF90 

Orifice  Ap  Transducer 

Transducer:  Statham  4347,  S/N  P6TC-2D-350,  ±  2  psig 
Signal  Conditioner:  Endevco  4423,  S/N  BA46 
Power  Supply;  Endevco  4425,  S/N  AF68 

Valve  Transducer 

Transducer:  Endevco  8510-B5,  S/N  PP75,  5  psig 
Signal  Conditioner:  Endevco  4423,  S/N  BA47 
Power  Supply:  Endevco  4425,  S/N  AF68 
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Xiklttt  (F8)  Transducer 

Transducer:  Endevco  8510B-5,  S/N  60HB,  5  psig 
Signal  Conditioner:  Endevco  4423,  S/N  BA44 
Power  Supply:  Endevco  4425,  S/N  AF70 

Exit  (PR4)  Transducer 

Transducer:  Endevco  8510-B5,  S/N  B34Y,  5  psig 

Signal  Conditioner:  Endevco  4423,  S/N  BA41 

Power  Supply:  Endevco  4425,  S/N  AF70 

Flow  Conditioning  Section  inlet  (PR2)  Transducer 

Transducer:  Endevco  8510-B5,  S/N  PP79,  5  psig 

Signal  Conditioner:  Endevco  4423,  S/N  BA63 

Power  Supply:  Endevco  442  ,  S/N  AF85 

Flow  Conditioning  Section  Exit  (PR3}  Transducer 

Transducer:  Endevco  8510-B5,  S/N  90EK,  5  psig 

Signal  Conditioner:  Endevco  4423,  S/N  BA63 

Power  Supply:  Endevco  4425,  3/N  AF85 

Section  1  Meem  Pressure  (Pmean)  Transducer 

Transducer:  Endevco  8510-B100,  S/N  80MB,  100  psig 

Signal  Conditioner:  Endevco  4423,  S/N  BA44 

Power  Supply:  Endevco  4425,  S/N  AF70 

Function  Generator:  Wavetek  142-S-527,  S/N  6700060 

House  Air  Temperature:  Thermocouple,  OMEGA  DP41-TC 

Atmospheric  Pressure:  Transamerica  Delaval  Digital  Barometer 

Type  2500-0103,  S/N  258 

Atmospheric  Temperature:  Bulb  thermometer 
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Frequency  Counter:  Racal-Dana,  Model  1992,  S/N  990308 
Valve  Mass  Flow  Meter:  Omega,  FMA-408-V,  S/N  7215 
Solenoid  Valve:  Automatic  Switch  Coirpany,  Model  8320A15,  S/N 
C78489S 
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The  errors  associated  with  the  dynamic  pressure 
measurements  are  dependent  on  the  accuracy  of  the  transducer, 
aitplifier  (signal  conditioner),  and  the  digitizing  card  in  the 
Nicolet  MultiPro  data  acquisition  system.  These  errors  are 
tabulated  below  from  information  provided  in  the 
manufacturer ' s  brochures . 

Transducer  (5  psig  Endevco)  Accuracy  =  ±  .2% 

Signal  conditioner  (Endevco)  Accuracy  =  ±  .5% 

Nicolet  MultiPro  digitizer  card  Accuracy  =  ±  45  uV 

For  gain  calculations,  the  pea)c  voltage  at  the  frequency 
of  interest  was  determined  from  the  spectral  analysis  program. 
A  representative  peak  voltage  from  the  gain  data  is  1.0  mV. 
For  the  1.0  mV  peak  signal,  the  actual  signal  from  the 
transducer  is  50  yV  since  the  signal  conditioner  amplifies  the 
signal  20  times.  Thus,  assuming  a  transducer  signal  of  50  yV, 
subject  to  the  accuracy  of  the  transducer,  the  resulting 
transducer  signal  will  be: 

(50  yV)  X  (±  .2%)  =  50  yV  ±  1  yV. 

When  amplified  by  20,  the  transducer  signal  and  error  become: 


(50  uv  ±  1  mV)  X  (20)  =  1000  mV  ±  20  uv. 


This  value  is  finally  subject  to  the  Nicolet  accuracy  of  ±  45 
UV,  yielding  a  final  range  of  actual  voltages  of: 

1000  uV  ±  65  uV. 

Thus,  a  1  mV  peak  displayed  in  the  spectral  analysis  program 
could  actually  be  as  small  as  935  uV  or  as  large  as  1065  uV, 
yielding  an  overall  error  in  the  dynamic  pressure  measurement 
of  ±  7%. 
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Several  programs  work  off  of  a  main  program  called 
maincrf.m.  The  main  program  and  all  of  the  subprograms  are 
included  in  the  following  pages.  A  listing  of  the  programs 
included  is  provided  below. 


maincrf.m .  E-2 

crf.m .  E-3 

exppipe.m . E-20 

FSpipe.m . E-25 

imp.m . E-27 

impa.m  . E-27 

adm.m . E-28 

adma.m  ...  E-28 

PENDPIN.m . E-29 

ZC.m . E-31 

dataA.m . E-32 

dataB.m . E-33 

dataBF.m . E-34 

dataBM.m . E-34 
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%  OMEGANUl,  OMEGANU2(J),  OMEGANU3,  OMEGANUB(B)  =  Characteristic 
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Units:  rad/s. 

•  for  valley (x,y)  array.  Units:  none, 
for  peak(x,y)  array.  Units:  none. 
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these  inputs . 

ZC.m  =  Accepts  input  variables  La,  Ca,  VELOCITY,  Z,  and  Y  and  outputs 

the  characteristic  impedance  for  the  individual  line  section 
as  adjusted  for  flow  effects. 
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INTERMEDIATE  VARIABLE 


AREA(I)  =  Area  of  the  inlet  (1),  flow  conditioning  section  (2),  and 
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function  Z  =  imp (OMEGARAT,  La,  OMEGA) 
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function  Y  =  ac3m{0MEGARAT,  Ca,  OMEGA,  NP,  GAMAIR) 
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function  thru  =  PENDPIN (VELOCITY,  La,  Ca,  ZC,  GAMMAL,  ZL) 
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Function  ZC.m  calculates  the  characteristic  iitpedance  for  a  line  with 
thrughflow.  This  function  is  called  by  program  crf.m. 
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thruzc  =  sqrt (Z/Y) * ( (sqrt (1  +  (m^2)))  -  m) ; 


set  A.  Full  scale  CRF  with  flow  straightening  elements. 
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%  Data  set  B.  Model  scale  CRF  with  flow  straightening  elements. 
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